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ABSTRACT 


Many problems in industrial quality control involve a vector of measurements of 
several characteristics simultaneously rather than a single characteristics For 
example, we may want to simultaneously control both the length and the diameter 
of the pipe In other words, both the length and diameter must be acceptable for the 
pipe to be usable Although one could monitor the piocess by separate control chart 
for each quality characteristics, but results of monitoring depend on the correlation 
between characteristics. If the conelation is high, the result may not be acceptable. 
So, controlling each characteristic separately may not yield a product m which all 
variables should be acceptable In many industries, separate control chart for each 
variable are used due to computational intricacies of multivariate control charts 
But, now with modern computers, one can easily use several types of multivariate 
charts. In this thesis, the peifonnance of different charts (some of which are 
multivariate charts and some are univariate charts) are computationally evaluated at 
different correlation coefficient, sample-size and process mean or process variance 
deviations levels Comparison is also done at different levels of three factors For 
different charts, different factors and their interaction are significant 
Pciformance of charts which aie basically for detecting process mean deviation or 
process variance deviation are also examined when process variance deviation or 
process mean deviation is introduced The present study thus provided some 
piehminary findings, which indicate that correlation coefficient, sample size and 
process mean or process variance deviation have effects on the performance of 
multivariate quality-control charts. 
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CHAPTER X 


INTRODUCTION 

There are many situations in which the simultaneous monitoring or control of two or 
more related quality chaiactenstics is necessaiy Foi example, suppose that a beaiing has 
both an innei diameter (x,) and an outei diameter (x 2 ) that togethei deteimine the 
usefulness of the pait Suppose that x, and x 2 have independent nonnal distributions As 
both quality chaiactenstics are measuiements, they could be monitored by applying the 
usual x chait to each chaiacleristic The piocess is consideied to be in control only if the 
sample means x, and x 2 fall within then lespective contiol limits 

Momtonng these two quality chaiactenstics independently can be very 
misleading The piobabihty that either x, 01 v 2 exceeds 3-sigma contiol limits is 

0 0027 However, the joint probability that both vaiiables exceed then control limits 
simultaneously when they aie both in control is (0.0027)(0 0027) = 0 00000729, which is 

considerably smaller than 0.0027 Furthermore, the probability that both x, and x 2 will 
simultaneously plot inside the control limits when process is really in control is 
(0 9973)(0.9973) = 0 99460729 Theiefore, the use of two independent x charts has 

distorted the simultaneous monitoring of x, and x, , in that the type I enor and the 
probability of a point eoirectly plotting in control are not equal to their advertised levels 
foi the individual control charts. 

This distortion m the process monitoring procedure increases as the 
number of quality characteristic increases. In general, if there are ^-statistically 
independent quality characteristic for a particulai pioduct and if X_Bar chait with P (type 

1 error) - a is maintained on each, then the true piobability of type 1 error for the joint 
control procedure is 

cc' =l-0-a) p 1.1 

and the probability that all p means will simultaneously plot inside their control limits 
when the process is in control is 


l 




P {all p means plot in control} = (1 - a) p 


1 2 

Clearly, the distortion in the joint contiol pioceduie can be severe, even foi moderate 
values of p Furthermore, if the p quality characteristics are not independent, which 
usually would be the case if they lelate to the same product, then equation-(l) and (1 2) 
do not hold, and we have no easy way even to measuie the distortion in the joint control 
procedure 

A Study on Multivariate Quality-Control Charts is done in the present woik In the 
next section pioblem definition, literature-ieview and the scope of present work aie 
discussed 


1.1 Problem Definition 

Process monitoung problems in which several related variables are of interest are called 
multivariate quality control problems. These processes can not be measured by using 
univariate control charts due to above-mentioned problem. Usually in the industlies the 
multivariate piocesses are handled by using univariate control charts due to 
computational and mathematical intricacies. But now with the help of modem 
computers, multivariate control charts can be used with ease 

In this thesis attempt is made to study different types of multivariate quality-control 
charts with iheii performance evaluation also. 


1.2 Literature Review 

The original work in multivariate quality was done by Hotelling (1947), who applied his 
procedure to bombsight data during World War II. He suggested the use of the statistics 



1 



A major disadvantage of Hotelling T2 statistics is that it does not indicates which variable 
or set of variables is out-of-contiol, when theie is signal 

Mason, Tracy, Young (1995) had piovided the approach foi identifying the souice of 
signals given by a T2 statistics in a multivanate contiol chait They decompose the T2 
statistics into independent paits, each of which is similai to an individual T2 variate, each 
of which reflects the contiibution of an individual vanable Each T2 is compared with 
some F distnbution value and based on this, decisions are taken 

Jackson (1980, 1991) discusses decomposing the T2 statistics into a sum of p principal 
components and using these components foi identifying the out_of_contiol variable 
He decomposed Z« as Y (l =C’AC' wheic C is the matiix of eigenvectors and A the 
diagonal matn\ of eigenvalues, the veetoi of scaled pnncipal components conesponding 
to ease X is given by U = A"^(’(A' -//„) Since T~ = U l U , so U is a decomposition of 
T2 into independent N (0,1) components 

When the pnncipal components lcpiesent meaningful gioupmgs of variables, the 
inteipretation of out_of_control signals is readily apparent However, in many cases it is 
difficult or even impossible, to attach meaning to the pnncipal components, and the 
characteristics associated with the oul_of_conliol signals can not be deteunined 

Jackson and Mudholkar (1979) also lecommended a procedure in which the leading 
pnncipal components of each case are computed and to predict the case 

Neduniaran and Pignuticllo .Jr. (1996) have also pioposed approach for diagnosing 
signals from 7'~ and control chaits They pioposed that once a signal is given on a 
T~ or a x 2 chart, univariate diagnostic charts are constructed for the p variables, any 
meaningful linear combination of the p vaiiables and then p principal components in 
such a way that the overall probability of false alarm is controlled at the specified value 
a for all the univariate control charts. When a signal is issued by a T 2 or a chart, 
statistics from the signaling subgroup are plotted on the 2p-univanate diagnostic charts to 



identify the variables and/or the pnncipal components that caused the signal Thus the T 2 
or a x 2 chart is used to monitor the multivanate process, whereas the univariate 
diagnostic charts are used only for diagnosing signals issued by the T 2 or a x 2 chart 

Doganaksoy, Tucker and Faltin (1991) suggested ranking the p variables based on their 
univariate t-statistics and find the P-value for each of the p variables when the null 
hypothesis of no change m their respective means is tested If the P-value of any variable 
is smaller than the corresponding Bonferrom error rate alp, one is to conclude that the 
mean of that variable has changed This process is similar to using p univariate control 

charts with limits x k ±t dj a/2p -\j s u ! n wheie x k and s kk are the standard mean and 
deviation of kth variable t dj ^ is the (1 - a I2p) th percentile point of a t-distribution 
with df degrees of fieedom, and n is the subgroup size 

According to Hayter and Tsui (1994), a “good” solution to the multivariate control 
problem is that the procedure should satisfy the following three conditions 

• Control the overall family wise enor rate at the nominal level a 

• Provide a simple mechanism for detecting which variables contributed to the 
out_of_control condition. 

• Quantify the magnitude of the mean change for the out_of_control variables 

They suggested the procedure besed on the construction of exact simultaneous 
confidence intervals for each of the variables means that meet each of these three goals. 
Both parametric and nonparametric piocedures are considered. 

In their approach, a signal is issued for the multivariate process when an M statistic 
exceeds the critical values C R>a The M statistics is defined by 


M = max 
!£/<£/? 


jy -A, 
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where X \ is the current observation of the ith variable, and /.gander , are the mean and 
the standard deviation of the ith variable, respectively When a signal is given, one has to 
examine each of the p univariate control charts with control limits set at //, ±C Ra cr l for 
diagnosing the signal 

Nedumaran and Pignatiello Jr. (1996)’ approach for diagnosing signal given by the T 2 
or a % 2 chart identified the variables or the principal components responsible foi the 
signals even m the cases where the Hayter and Tsui (1994), method failed to signal 

Holmes and Mergen (1993) have also pioposed technique for out-of-control signal on 
the T 2 or % 2 chart 

They proposed to construct standardized value for each of the p variables when there is 
signal The standardized values are then compared with control limits of 0 ±Zy 

The limits are calculated using the Bonferrom appioach which says that the effective 
value of a (l e , Type I error) risk (eff a) is approximately 
eff a =pa , 

Where a is the risk associated with each variable 

This method is different fiom the other methods in that capability standard deviation is 
estimated using the mean squaie successive difference (MSSD) Also standardized value 
are plotted simultaneously on the same display 

Timm, N.H. (1996) suggested Finite Intersection Test (FIT), m which they use single 
step and stepdown FIT to evaluate whether a multivariate process is m control or 
out_of_control The FIT computer piogram is also developed by Cox, Frag, Boudreau, 
and Timm (1994) . 

The stepdown FIT procedure provides an alternative to the procedure proposed by the 
Hayter and Tsui (1994), the sequential test procedure by the Jackson and Bradley 
(1996), and regression adjustment procedure of Hawkins (1991, 1993). The procedure is 
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more powerful than Roy’s stepclown recommended by Wierda (1994b) and the 
procedure of Mason, Tracy, and Young (1995) 

When multivariate normality does not hold, the stepwise resampling method of multiple 
hypothesis listing developed by Troendle (1995) may be used. Troendle’s procedure 
incorpoiates the stiucture of the covariance matnx into the problem without making 
distribution assumptions Howevei, it does not provide simultaneous confidence interval 
for population paiameteis 

In multiple processes, a Shewhart chait is often used at each individual step If the steps 
of processes are independent then this is a meaningful procedure However, when steps 
aie correlated, it is difficult to interpret Hotelling T 2 can be used here, but process 
quality charactenstics should be normal landom variables Also, when the process is out- 
of-control, it is difficult to say, which variable is responsible 

Zhang (1980, 1982, 1984, 1985a, 1985b, 1989a, 1989b, 1992) had proposed a different 
approach foi this Zhang (1984) call his chart “cause-selecting chart”, which is similai to 
in concept to the regression contiol chart of Mandal (1969) m that a control chart is 
constructed for a vanable only after the observation have been adjusted for the effects of 
some other landom variable Cause-selecting charts use incoming quality measurements 
ant out-going quality measurements m an attempt to distinguish between m-coming 
quality problems and problems in the current operation It has advantage ovei the 
multivariate T 1 chait in that it provides diagnostic information regarding which 
subprocess is out of control Also it handles complicated relationship between the 
variables. 

Wade and Woodall (1993) propose using prediction limits with cause-selecting chart to 
improve their statistical performance 

When false-alarm rates of the cause-selecting chart using prediction limits and Zhang’s 
limits were examined, it was found that the false-alarm rate of Zhang’s limits depends on 
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the sample size foe estimating the linear relationship Thus the advantage of using 
prediction limits is that the false-alarm rate is not affected by the sample-size 

Anothei chart using T 2 is that of Alwan (1986), who proposed a cumulative sum 
(CUSUM) of the T 2 

Crosier (1988) suggested a CUSUM of T (COT) But CUSUN of T2 is preferred over 
CUSUM of T due to computational simplicity 

Crosier also proposed another CUSUM scheme, refereed to as multivariate CUSUM, 
forms a CUSUM vector directly from the observations CUSUM vector is not based on 
the T 2 of the individual observation, but its signals are based on a Hotelling T 2 applied 
to the CUSUM vectois 

Both procedures allow use of recent enhancements for the CUSUM schemes, such as FIR 
feature Both of this CUSUM procedure leduces to multivariate Shewhart charts when 
the CUSUM scheme parameter h is 0 

Another approach to multivariate CUSUM control is the proposal of Woodall and 
Ncube (1985). This consists of a set of individual CUSUM’s of the variables of X, and a 
signal is given if any of the p CUSUM’s exceeds its control limit This approach has the 
bettei mterpretabihty, since signal obtained corresponds to an immediately identified 
variable But, when vanables are highly correlated, then this gives poor result 
In the case of strongly correlated variables, they suggested the use of principal 
components for improvement of the performance Though signal is given for a 
component rather than a particular variable, there may be some loss of mterpretability m 
term of the original variables. 

Pignatiello and Runger (1990) also pioposed two multivariate CUSUM charts 
First multivariate CUSUM scheme is based on the following vectors of cumulative sums 


A - ±x. 

/=/-/!+! 
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and 


MC t = max jo, (d'Z x D l Y -& 2 / ( 

where k 2 > 0 and 

/,=/,_,+ 1 if MC,_ t > 0 

=: 1, otherwise, 

I = 1,2 ,3, an out-of-control signal is given as sonn as MC l > h 2 where h 2 >0 

Second CUSUM chart is the CUSUM of Hotelling T 2 

First multivariate CUSUM scheme of Crosier is similar to this, but yet distinctively 
different from this CUSUM vector scheme 

Crosier “contracts” or “shrinks” each of his vectors of cumulative sums toward the zero 
vector by multiplying the cumulative sum by a scalar The norm of the contracted 
cumulative sum is then compared with an upper control limit The ART performance of 
Crosier’ multi vauate cumulative sum chart is similar to that of Pignatiello and 
Runger’s. 

Woodall and Ncube’ multiple univariate CUSUM scheme and Pignatiello and 
Runger’s CUSUM vector scheme give roughly the same performance when the process 
experience a mean shift in just one variable or characteristic But when means of several 
variables shifts simultaneously to off-target values, univariate CUSUM chart gives poorer 
result then CUSUM vector scheme 

Lowry, Woodall, Rigdon (1992) proposed multivariate extension of the univariate 
exponentially weighted moving average (EWMA) chart The ARL performance shows 
that MEWMA procedure can perform better than the multivariate CUSUM procedure of 
Crosier (1988) and Pignatiello and Runger (1990) when the process is initially out-of- 
control and it performs roughly the same if the shift in the mean vector is delayed. 
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Healy (1987) provided a geneial mullivanate framework to produce more powerful 
diagnostics than those based on the T 2 statistics 

Suppose it is anticipated that when the process mean shifts from to jU A Then the 
optimal test statistic for a shift m the ith measured vector is the scalar quantity 

(Ma-M)' r'(X,- J u) = (u A - J u)'Z l 

wheie the p x 1 vector Z, is defined by 

Z, =Z' , (X ; -;0 

wheic vector Z ; is actually a vector of regression-adjusted variables 

The ith component of the Z ( is the standardized residual resulting when X, is regressed 

on all other variables in X 

Hawkins (1991) put Healy’s result and Mandal (1969) regression adjustment together 
and pioposed a different appioach 


1.3 Scope of Present Work 

The present work focuses on the following issues 

• To study the effect of correlation coefficient between different variables, sample size 
and deviation in the process mean or process variance on the different charts 

• To study the effect of mean deviation on the charts meant for measuring process 
variance deviation and the effect of variance deviation on the charts meant for 
measuring process mean deviation 
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1.4 Organization of The Thesis 


Chapter 2 entitled ‘Description of Multivariate Quality-Control Charts’ discusses 
different control charts for measuring process mean and variance deviation It also 
discusses the merits and demerits of these charts 

Chaptei 3 entitled ‘Performance Evaluation’ covers two issues namely performance 
of different charts with respect to con elation coefficient between different variables, 
sample size and deviation in the process mean or process variance on the different charts 
and the comparison of different charts 

Chapter 4 entitled ‘Results and Conclusions’ summarizes all the findings that have 
been made through out the work and the conclusion that follow from these findings 
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CHAPTER 2 


DESCRIPTION OF MULTIVARIATE QUALITY-CONTROL 
CHARTS 


This chaptei discusses different types of quality control charts Some charts are only for 
momtonng process mean changes and some are only for monitoring variance changes 
Both types of charts are discussed Their merits and demerits are also focused 

Charts for Monitoring Mean Changes in The Multivariate Processes: 


• CHI-SQUARE ( x 2 ) CONTROL CHART 

• T SQUARE (T 2 ) CONTROL CHART 

• MULTIVARIATE CUSUM CHARTS : 

(1) CUMULATIVE SUM OF T 

(2) CUSUM VECTOR SCHEME 

• MULTIVARIATE EXPONENTIALLY MOVING-AVERAGE CHART 

• REGRESSION-ADJUSTMENT PROCEDURE 

(1) Z-METHOD FOR MONITORING PROCESS MEAN 

(2) Y-METHOD 

• UNIVARIATE SHEWHART X-BAR 

• UNIVARIATE X-CUSUM CHART FOR MONITORING PROCESS MEAN 


Charts for Monitoring Variance Changes in the Multivariate Processes: 


• REGRESSION-ADJUSTMENT PROCEDURE 

1. Z-METHOD FOR MONITORING PROCESS VARIANCE 

• UNIVARIATE X RANGE CHART 

. MULTIVARIATE PROCESS VARIANCE CHART 

• UNIVARIATE S 2 CONTROL CHART 



• UNIVARIATE X-CUSUM CHART FOR MONITORING PROCESS 
VARIANCE 

2.0 MONITORING MEAN CHANGES IN MULTIVARIATE 
PROCESSES : 

2.1 CHI-SQUARE ( Z 2 ) CONTROL CHART 
Assumption: 

The joint probability distribution of the p quality charactensrtic is the p variate normal 
distubution 

This procedure requires computing the samples mean for each of the p quality 
characteristics from a sample of size n This set of quality characteristic mean is 
represented by px 1 vector 



Where a is the confidence level 

The test statistics plotted on the Chi-Square Chart for each sample is 

j 2 =n(x- / i/)'E _1 (x-^) 211 

Where p =\ju x ,ju 2 , . .,p p \ is the vector of m-control means for each quality 
characteristics and E is the co-variance matrix The upper limit on the control chart is 

UCL = 2 1 2 
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Estimating fj. and E 


In practice, it is usually necessary to estimate // and E from the analysis of preliminary 
sample of size n, taken when the process is assumed to be m control If m such samples 
are available the sample means and variances are calculated from each sample as 
follows 




1 V- 

fj = 1=2, 

»p 

2 1 3 

x Jl = 

~Z x ur 

[k =1,2, 

>p 


s 2 - 

_ i !>,, x ir) 2 | 

a i 1=1 | 

fj = U. 

[k = 1,2, 

»p 

,p 

2 1 4 


Where x k is the ith observation on the jth quality characteristics in the kth sample. The 

covariance between quality characteristic j and quality characteristics h in the kth. sample 
is 


s ihk -XjkXXhk - x hk) 

n -1 1=1 


Jk =1,2, ,m 


2.1.5 


The statistics x jk ,S 2 ]k and S jhk are then averaged all over m samples to obtain 


and 


= 1 ^ 
x i = ~2, 

m tt 

i m 

s?=-2s 

mti 


k jJv 


1 m 

S *=-E S * 

mt; 


j^h 


J 1j 2 3 . . ,p 
j = 1 >2, ,p 


2.1 6 

2.1 7 

2 1 8 
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The Xj are the elements of the vectoi X, and the p x p averages of sample covariance 
matrices S is formed as 


S = 






2 1 9 


The averages of the sample covanance matrices S is an unbiased estimate of 2 when the 
process is in control 


2.2 T-SQUARE (T 2 ) CONTROL CHART: 

Assumption: 

When the process is in control, the X, are independent and follow a multivariate normal 
distribution with mean vector p and covariance matrix 2 


If we replace p with x 


and 2 with S in Eq-2 1 1, test statistics now becomes 


T 2 



22 1 


In this form , the procedure is called the Hotelling T 2 Control Chart 


The control limit for T 2 control chart are given by 

UCL = 

LCL = 


p(m-l)(n-l) 

p, mii-in-p + 1 

mn - n - p + 1 
0 


2.2.2 
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2.2.1 The Case Where n = 1 


In some settings the subgioup size is generally n = 1 
Example : chemical industries 

Suppose that m samples, each of size n=l are available, and that p is the number of 
quality characteristics observed in each sample Let x and S be the sample mean vector 

and covariance matrix, respectively, of these observations The Hotelling T 2 statistics in 
Eq-2 2 1 becomes 

T 2 = (x-x)'S"'(x-x) 2 2.3 


The control limits are given by 

Ucl = P( m + l)( m ~l) r 
m 2 - mp 

LCL = 0 


224 


when the number of sample m is large , say m > 100 , upper control limits is given by 


or 


UCL = 


P(m-l) 
m - p 


a,p,m-p 


UCL = xl. P 


225 


226 


2.2.2 Merits and Demerits T 2 or j 2 Control charts 

Hotelling T 2 and % 2 control charts are two of the most popular multivariate control 
chart, because of its similarity to a univariate Shewhart control chart It is often called 
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multivariate Shewhart- type control chart These charts are widely used when the 
observations are independent and identically distributed p-variate normal random 
variables The probability of incorrectly declaring the process to out of control can be 
maintained at any desired level in these charts 

These aie the optimal test statistics for detecting a general shift in the process mean 
vector for an individual multivanate observation, but it is not optimal for more structured 
mean shifts — for example when shifts m only some of the variables 
The major drawback of this chait is that they are unable to predict which variable of set 
of variables is out-of-conti ol Also mean shifts are always confounded with variance 
shifts, which require quite extensive analysis following a signal to determine the nature of 
the shifts These proceduies are based on only the most recent observation, it is 
insensitive to small and moderate shifts in the mean vector 


2.2.3 Identifying The Sources For Out-of-Control Signals When The T 2 or 
X 2 Control Chart is Used: 

Difficulties encountered with any of the multivariate control chart are practical 
interpretation of an out-of-conti ol signal Which of the variable of (set of variables) is 

lesponsible for signal is difficult to say The standard practice is to plot univariate x 
charts on the individual variables x, , x 2 , . ,x p But this approach may not be 

successful due to correlation of variables 

T 2 or x 2 value which is out-of-control indicates that one or more of the mean has 
changed or that the correlation between two or more of the variables has changed, but its 
difficult to determine which valuable is out-of-control 

Following method is one of the recommended methods for identifying the out-of-control 
variable 

Donald S. Holem and A. Erhan Mergen (1993) had used Bonferrom-type control 
limits The method is done in two steps as follows 
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In first step, we calculate the standardized value for each of the p variables The 
standardization is done like this 


-r£ 


m ~ 


22 3 1 


Wheie m is the number of samples s t is the mean square successive difference 
estimatoi for the capability standard deviation It is defined as 


fig 

zA+i ~ x ')‘ 


= 


(=1 


2(u-l) 


( 3 

1 - — 

8 n 


v-i 
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In second step, we make a control chart type of display of the standardized values 
The limits on the display are calculated using the Bonferrom approach which says that 
the effective value of the a (i e , type I error) risk (effective a ) is approximately 

eff a = pa 2.2 3.3 

where a is the risk associated with each variable 

The approximate z for the control limits is established using a = eff a /p Then the 
limits are set at 0 ± z a/ 

n 

Example: 

When p = 2 and eff a = 0 003, the corresponding a and z values 0 0015 and the 
z omn5/ , =3 19 respectively The variable plotting outside this control limit is considered 

responsible for signal If all standardized values fall within control limits, the correlation 
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of variables may changed This method is different fiom the othei methods m following 
aspects The first difference is that we estimate the capability standaid deviation using the 
mean square successive difference MSSD approach The second difference is the 
simultaneous plotting of the standardized values on the same display 


2.3 MULTIVARIATE CUSUM CHARTS : 

The most frequent application of the CUSUM schemes is the detection of change m the 
mean of a normally distributed variable 

Ronald B. Crosier (1988) had given two-multivariate cumulative sum (CUSUM) quality 
procedure The first CUSUM pioceduie reduces each multivariate observation to a scalar 
and then forms a CUSUM of scalais The second CUSUM procedure forms a CUSUM 
vector directly from the observations 


2.3.1 CUSUM # 1 : 

The most direct method is to form a CUSUM of scales T n (n = 1,2,3, ) This method is 

called COT scheme 
The CUSUM is given by 

S„ = max(0,xS' ( „_ 1) +T u -k) 2 3 11 

where S 0 > 0 and k>0 

T n ( Square root of Hotelling T 2 ) is given by 

2 3.12 
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The COT scheme signals when S n > h (Upper Control Limit) 

The COT scheme offers several advantages over multivariate Shewhart charts First COT 
scheme may be designed to give any false alarm rates by using desired type I error 
Second, recent enhancements for CUSUM scheme may be applied to COT schemes 
A lobustness procedure, such as the two-in-a-row rule (Lucas and Crosier (1982)) may 
be applied to the sequence of T statistics The two-in-a-row rule states that outliers should 
be ignored (1 e not used m the CUSUM calculations) but the occurrence of two outliers 
in succession is taken as an off-aim signal The fast initial response (FIR) feature for 
CUSUM schemes (Lucas and Crosier (1982)) is easily extended to COT schemes 
Quick detection of an initial off-aim condition is obtained by starting the scheme with S Q 
equal to h/2 lather than 0 


2.3.2 CUSUM #2 : 

A vector valued CUSUM scheme is derived by replacing the scalar quantities of a 
univariate CUSUM scheme by vectors 
The univariate CUSUM scheme is given by 

S„ = max(0, S {n _ X) + (X„ - pi)-ko ) 2 3 2 1 

where p is the taiget value for the mean, a is the standard deviation of the X’s, k>0 and 
S 0 =0 

Replacing the scalars by vectors gives s„ = max(0, s n _ { + (x - pi) - k ) . 

But the problem is then, how to find k and how to interpret taking the maximum of vectoi 
and the null vector In the univariate case, the quantity S (n _ {) + (X n - pi) is shr unk toward 

0 by k standard deviation But for multivariate case k must satisfy k'Y,~ [ k = k 2 , that is k 
must be of length k, where length is defined using the covariance matrix 2 If the 
subtraction of k is to shrink to +x n ~ pi toward 0, then k must be in the same 
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direction as s n _ x +x n -/j. Hence k = (k/ C„ + x n - pi ) , where C n is the length 

+ x n - j-i The maximum of null vector can be taken as s„ = 0 , rather than 
subti acting k, whenever the length of k is greater than the length of s n _ x + x n - ju, that is 
whenever k>Cn Rather than calculating the vector k and subtract it from +x n - ju 
It is simpler to just contract the vector + x n - ju by (1-k/Cn), provided that k<Cn 
Hence multivariate vector CUSUM method is expressed as 

C n = [(Vl + X n - /0'S" 1 (Vl + X « - /“)] /2 2 3 2 2 


then 


s 


n 


= 0 


if C„<k 23 23 


= ( V. +x„-M -k/C„) If c n >k 23 24 


where s Q = 0 and k>0 



The multivariate CUSUM vector scheme signals when Y n > h 


23 2 5 


The multivariate CUSUM vector scheme give quicker detection of small shifts than 
both multivariate Shewhart charts and COT schemes when n = 1. But, for n>l, 
COT method gives better result than Cusum vector method. The multivariate 
CUSUM chart also don’t give the indication that which variable is out-of-control 
when there is signal 


2.4 MULTIVARIATE EXPONENTIALLY WEIGTHED MOVING-AVERAGE 
CONTROL CHART : 
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The exponentially weighted moving-average (or EWMA) control chart is also a good 
alternative to the Shewhart chart when we are interested in detecting small shifts The 
perfoimance of the EWMA control chart is approximately equivalent to that of CUSUM 
chart, and in some ways easy to set up and operate It is also commonly used with sample 
size n=T 

The univariate EWMA chart is based on the values 

Z, =rX, + (l-r)Z,_, 24 1 

i=l,2, where Z 0 = ju 0 = 0 and 0 < r < 1 

If Xi are independent N(0,<r 2 ) random variables, then the mean of Z, is 0 and the 
variance is 

a\ t - {i'[l - (1 - r) 2, ]/(2 — r)]cr 2 , 1 = 1,2, 

Thus when the in control value of the mean is 0, the control limits of the EWMA chart 
are often set at ± Z.cr z , where L and r are the parameters of the chart 

In the multivariate case, a natuial extension is to define vectors of EWMA’s 

Z, =RX, +(/-i?)Z M 242 

I = 1,2,3, , where Z 0 = 0 and R = diag( r ]} r 2 ,. r p ), 0 < r < 1, j=l,2,3, ,p 

The MEWMA chart gives out-of-control signal when 

T, 1 = ZX Z 'Z, > h 243 

where h (>0) is chosen to achieve a specified incontrol ARL and Z Z; is the covariance 
matrix of Z, 

Z z is given by 
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244 


£z, = 't[R(I-R)'- J Al-Rr , R\ 

/= i 

if r, = r 2 = = r p = r, then 


2 Z is given by 


2 Z , ={r[l-(l-r) 2 ']/(2-r)}s 2 4.5 

if there is no a prion reason to weight past observations differently for the p quality 
characteristics being monitored, then r, = r 2 = =r p = r 

if r { = r 2 = = r p = r, then the MEWM A vectors can be written as 

Z, =rX, +(l-r)Z,_, 24 6 

Smaller values of r are more effective in detecting small shifts m the mean vector 

The ARL performance of MEWMA chart is similar to_that_of multivariate CUSUM 
chart. But when n > 1, we find earlier signal in multivariate CUSUM chart than 
MEWMA. This chart also doesn’t have provision for identification of out-of-control 
variables 

2.5 REGRESSION ADJUSTMENT FOR VARIABLES IN MULTIVARAIATE 
QUALITY CONTROL : 

Problems m multivariate quality problems occur when one observes, not a single 
characteristics, but a vector of characteristics It is of course possible to monitor the 
individual characteristics by charting them separately, but to extent that they are 
correlated. It can be more effective to monitor a vector of characteristics m some 
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multivariate way The mam point of i egression adjustment is that the nature of this 
multivariate focus should, wheie possible, take into account the characteristics 
mechanisms that lead to the association between the process measures Two important 
possible mechanisms are the following 

• A mechanism that leads to some or all of the measurement being related, but where 
possible foi any of them to undergo parameter shifts without affecting the others 

Example: 

measurements of different dimensions of switch drum it is possible here that 
one dimension can shift , without affecting otheis But response to problem can be 
speeded up by appiopnate use of the measures on the other dimensions 

• A mechanisms that leads to measures having a natural ordering in which if any 
measures undergoes a parameter shift it may affect some or all of the measures 
following it, but none of the measures preceding it m this ordering This sort of 
ordering is called a cascade And the resulting vector of measurements described as 
having the cascade piopeity 

Cascade property is typified by serial value-added manufacturing processes 

Example: 

A metallurgical flotation plant wheie the material passes through successive enrichment 
steps The quality of product at the end of any steps depends not only on that step is 
opeiatmg correctly, but also on how well the preceding steps performed If a particular 
step goes out of control, it can affect product measures, not just at that step, but also some 
or all of the following steps 

Only the first mechanism is discussed in this thesis. 

For the first mechanism, Hawkins (1991) defines a set of regression - adjusted variables 
using the vector 
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( _I YT* 

diagl^ Y = A(X-jU 0 ) 
v o )_ 

where the transformation matrix A is given by 

( -\ \ -i 

A = dictgl £ 

V o J J o 

Z, is the standaidized residual when the ith variable is regressed on the remaining p-1 
vauables in X 

In the second type of mechanism, Hawkins decomposes the T 2 statistic using the 
standaidized residuals from the regression of the ith variable on the first i-l variables 
These are defined by 

Y = C{X - ju) 253 

wheie C is the Cholesky lower triangular root of S' 1 

2.5.1 Improvement in Performance Due to Regression Adjustment 

The diagnostics Y and Z have two motivations One is that they provide controls on the 
scale of the individual measures and the other is that they are more sensitive to out-of- 
control situations because the regression adjustment reduces the background vaiiabihty, 
letting the shifts stand out more clearly The relative performance of regression-adjusted 
and unadjusted diagnostics depends on the degree of correlation on which the regression 
is based Suppose that on the original scale of the mean of X ] shifts by 6 standard 

deviations and that an adjustment is made by regressing X , on the variable that affect it 
Let this regression have a multiple conelation coefficient of R / The regression-adjusted 
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vanable will also have its mean shifted by 8 standard deviation of X / , but since the 
1 egression residuals have a standaid deviation smaller by a factor of 


This shift will be 
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of the regression-adjusted variable’s standard deviations At one extremepf Rj is 
zeio,then there is no improvement in performance over working on the x scale At the 
other extreme ,if Rj is close to l,the margin of improvement of the regies sion-adjusted 
vanable over the original variable can be arbitrarily large 

This is also important in regression adjustment is that, the variable for which regression- 
adjustment is made, they did not undergo any shift in mean m mean If they did, then the 
calculated improvement in performance would not occur This has implications for the 
type of regression adjustments made Suppose that measurements have cascade property 
Then it is correct to adjust each X f for A,, , X jA but not for any of 


X lA , , X p (i e one should use the diagnostics Y t ) Suppose Z ; are erroneously used 
instead, so that X is regiession-adj usted not only for X x , ,X jA , but also for 
X i+] , ,X , which are they affected by the shift m X } Then the regression- 
adjustment would incorrectly attenuate the mean shifts so that performance of Z j would 


not be as good as that using Y / and perhaps not even as good as that using X } itself 


2.5.2 REGRESSION ADJUSTMENT METHOD FOR TYPE 1 MECHANISM : 

Assumptions: 

• When the process is m control, the X ~ N(p,2) are independent and 
follow a multivariate normal distribution with mean vectoi p and 
covariance matrix E 
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While the piocess is in control, p = /u 0 and 2 = 2 0 


• All the components of X are rescaled to unit standard deviation — that is, 
that the diagonal elements of 2 0 are equal to 1 


Suppose <7 n and a 11 and are the i,jth elements of 2 0 and 2 0 1 respectively 
The multiple regression of X, , the ith component of X, on all other components of X is 
wntten as 

X, - H = 2 (X j - M t ) + £ , 2 5 2 1 


r„ = residual variance i e the variance of the e, 


In terms of variance and covariance 


T „ =&„ 

cr" =-J3 Jr n {i*j) 


2 5 2.2 


Instead of completely unspecified shift in p, assume that the out-of-control affects the 
mean of a single variable l, shifting the process distribution to N(ju 0 +Se t , 2 0 ), where 
e t is a column vector whose value is 1 in row 1 and 0 elsewhere and the scalar 8 is the 
magnitude of the shift then the likhhood-ratio statistic for the hypothesis 8 = 0 is given by 


Z, = 







2 5 2 3 
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Whose null distribution is N(0,1) This (scalar) Z, is the residual when X, 
is i egression all othei components of X lescaled to unit variance 

The vector of scaled residual is also expressed as 

Y,=±iX-Mo) 25 24 

0 

From Equations-2 5 2 4, the ith component of Y is the regression residual when variable I 
is legressed on all other variables, scaled by a factor t~ 1 . Under control, 7~ N (0, ), 

Z is just rescaling of Y 


Z = 


diag 


[ 1 ) 


1-Z 


Y = A(X -» 0 ) 


2 5 2.5 


where the transformation matrix 


A = 




J ° 


252 6 


The coefficients of A provide an instant visual diagnostics of the relative effectiveness of 
multivariate control using Z over univariate control using X If A is close to an identity 
matrix, then theie will be little magnification of mean shifts in going from X to z, so the 
proposed methods will not be noticeably more effective even for single-variable shifts 
than univariate control methods Against this, numerically large elements m A 
(particularly its diagonals) imply laige magnifications of mean shifts and large potential 
benefits from using this method The matrix A also shows the effects of changes m more 
than one elements of ju Q A vector shifts of say Ap in the mean would translate into a 

vector of shifts in Z of AAp, which is easily evaluated 
It follows from definition that Z ~ N (0,B), where 
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r > 

- 


( -1 \ 

B = 

N° 

So ^ 

Sj 

1 


Z 

0 

J 


The relationship between X, Z and T 2 is 



= Y j (X- j u 0 )Y i 2 52 8 

o 


oi, rewriting m terms of Z, 


T>=± W, 

/ = 1 

Where W t = (X, - /.i, )Z l r' A 2.5 2 9 


The vector Z provides the possibility of separate control chart of the individual variables 
in X , thus gives trigger for the out-of-control variable also 

The major advantage of Z_method is that we don’t need separate diagnostic procedure to 
know about the responsible variables for out-of-control sample But when we run the 
process even after out-of-control signal, the signals get mixed with other variable also, 
and undeviated variable also starts to give out-of-control signal Since one is usually 
interested both in special causes affecting the particular case and possible persistent shifts 
in the process, it is generally desirable to use combined Shewhart/CUSUM control charts 
In this thesis Shewhart X-bar, X-Range chart, X-Cusum for mean and scale for individual 
X’s and Z-Cusum chart for mean and scale for Z elements have been used 
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2.5.2.1 Cusum Control Chart for Z elements : 


CUSUM control of location and scale is given by as follows." 

Vo = Vo = S, + o - S- = 0 
W m =( ]Z,/-0 822)/0 349 

K = max(0, + Z,„ -k) 2 5 2 1 1 

Vn= mm(0, +Z m +k) 

Sl= max(0, 

=mm(0, 5“,., +lf m + *) 

The decision -interval CUSUM’ s Z + and L test for shifts of location m the upward 
and downward direction, respectively The CUSUM’ s S + and S~ test for shifts of 
variability m the upward and downward directions, respectively Their basis is the 

observation that if Z ~ N(0,1), then |Z|^ closely approximates an N (0 822,0 349 x 

0 349) distribution , and that changes in the scale of Z affect the location of |Z|^ By 

plotting all four CUSUM’s on a common chart, one obtains a powerful CUSUM control 
foi location and scale An out-of-control is signaled if any of these four CUSUM 
quantities exceeds the decision interval h Along with this Cusum chart of mean and 
scale, Shewhart X_bar chart and X_Range chart are used for X’s, which are described 
below 


2. 5.2.2 UNIVARIATE SHEWHART X-BAR AND X-RANGE CHARTS : 

Assumption: 

Quality characteristics are normally distributed with mean jU and standard deviation 
a These charts are used in Z_method of regression adjustment procedure This chart is 
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developed by Dr. Walter A. Shewhart (1920). We usually do not know fj. and a , so 
we estimate them from past data, which are m control But, when we have standard value 
of /u and a , we used these charts with standard value of fi and a 


Charts Based on Estimated Value of /.i and <x : 


X_bar chart: 

If x, ,x 2 , , x„ is a sample of size n, then the averages of this sample is 

- X, + X, + + x„ 

X = — - = - 


2 522 1 


n 


If x, , x, , x 3 , , x m are the averages of sample m, then best estimator of ju is given by 


x = ■ 


X, + x 9 + 


+ x„ 


m 


2.5 2 2.2 


Control Limits for X_bar chart 

UCL = x + A 2 R 

Center Line = x 2 5 2.2 3 

LCL = x-A 2 R 

Where A 2 is constant depends on sample size 


Range Chart: 

In this chart, we estimate standard deviation from ranges of the m 
samples The range of the sample is the difference between the largest and smallest 
observation m sample size n, that is 

^ Tna\ Trim 

If R ^ , R 2 , . , R m be the ranges of the m samples. The average range is 
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m 


2 5 224 


- _ R x + R 2 + +R n 


Control Limits for Range chart: 

UCL = D 4 R 
Centei Line = R 
LCL = D 3 R 

The constants Z) 3 and D 4 depends on sample size n 
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Charts Based on Standard Values of // And a : 

When we know the standaid vales of /./ and u , we use these standard to establish the 

control charts for x and R without analysis of past data 

Suppose that the standaids for mean and standard deviation is // and cr 


Control Limits for X bar chart: 


UCL = ju + 

■\Jn 


Center line = ju 


LCL 



Control Limits for Range chart: 

UCL = Oo- 
center Line = d 2 cr 
LCL = D x a 

Where constants D 2 ,d 2 , D ] depend on sample size n 


25 226 


2 5 2.2 7 
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If the sample size is relatively small, the range method gives almost as good an estimator 
of the variance a 2 as the sample vanance S 2 For modeiate values of n, say n > 10, the 
range method loses efficiency rapidly, as it ignores all the information m the sample 
between x mav and x min So we use S control chart for sample size > 1 

Control Charts for XBar and S: 

x and S aie used foi all values of n as it gives good result then range method Earlier 
range method was used due to simplicity of calculations In these charts we calculate 

sample average x and the sample standard deviation S for each sample If a 2 is the 
unknown variance of a probability distribution, then an unbiased estimator of o 2 is the 
sample variance 

£(*, -*) 2 

S 2 =-^ 

n - 1 

But the sample standaid deviation S is not an unbiased estimator of cr Actually, S 
estimates c 4 a , where constant c 4 depends on sample size n 

Charts Based on Standard values of /u and a : 

When we have standard value of /u and a , we don’t estimate it from past data 

Control Limits for the x chart: 

The parameters of x chart is given by Eq-2 5 2 2 6 

Control Limits for the S Chart: 

UCL = B 6 <t 

Center line = c 4 cr 2.5 2 2 8 

LCL = B 5 a 

Where constants B 6 ,c 4 ,B 5 depend on sample size n. 
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Charts Based on Estimated values of ju and a : 


If no standaid is given for x and a , we estimate it from past data Let S l is the 
standard deviation of the ith sample The average of the m preliminary samples is 


_ 1 m 

s-rZ®. 


m 


(=i 


Standard deviation a is estimated as 


Control Limits for the x chart: 


UCL = x -t- 


_3S_ 

c A Jn 


Centerline = x 


LCL 


x — 


3S 

c.-Jn 


25229 


Control Limits for the S Chart : 

UCL = B 4 S 

Center line = S 252210 

LCL = B 3 S 

Where constants B A and 5, depend on sample size n. 
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2.5.2.2.1 Control Charts for Individual Measurements : 


When the sample size n =1, we do not use above mentioned x and S charts or x and 
Range charts Here, we use moving range concepts 

The moving lange MR, is given as |x, - x,_, | , where x, is the observation of ith sample 
We establish conti ol charts for monitoring x and cr based on moving range 


Control Limits for the x Chart: 

UCL 


X + J - 


MR 


d. 


Center line = x 


LCL 


= x-3 


MR 


2 5 22 1 1 


Constants d 2 is taken as value conespondmg to sample size n - 2, since we substract 
consecutive observation 

Here, we use moving range chart for monitoring variance deviation 


Control Limits for Moving Range Charts: 

UCL = D 4 MR 

Center line = MR 2 5 2 2 1 2 

LCL = D Jar 

Heie, also we take value of D 4 and D-, corresponding to n = 2 

2.S.2.3 UNIVARIATE X-CUSUM CHART : 

This chart is also used in Z_method of regression adjustment procedure The major 
disadvantage of any Shewhart control chart is that it only uses the information about the 
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piocess contend in the last plotted point and it ignoies any information given by the entire 
sequence of points Due to this reason, the Shewhait control chart is sensitive to small 
shifts in the piocess For this reason we use cumulative sum control chart of small sets 
Because it considers all the information in the sequence of sample values by plotting the 
cumulative sums of the deviation of the sample values from a target value 
Cumulative-sums control charts weie first proposed by Page (1954). 

Assumption: 

The quality characteristics are normally distributed with mean // and standard 
deviation c 


If x, ,x 2 , , x m are the individual observation of m samples, the standardized two-sided 

Cusum is given by 

C, + = max 0, y, - k + C ( t, ] 

C~ = max 0 ,-k -y, + C~_ { 



where, 

If sample size n >1, 


y, = 


*, -Mo 

<T 


y, = 


x , -Mo 
a I yfn 


and Cq = Cq = 0 

Where, x, , x 2 , , x m are the averages of m samples of size n and and C 0 + = C 0 ~ = 0 

k = K/cr 

K is called the reference value and is given by, 


where, //, = y/ 0 + Scr = deviated mean from target value 

If either C, + or C~ exceed the decision interval h, the process is called out of control 
Due to standardization many Cusum charts can use the same value of k and h and these 
parameters do not depend on a . 
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2 5.2.3. 1 Cusum for Monitoring Process Variability: 


Cusum chart can also be used for monitoring process variability 


Hawkins (1981)(1993a) has suggested a new standardized quantity 

^|-0 822 


0 349 


2.5 2 3 1 1 


He suggests that the v, are sensitive to various changes rather than mean changes But m 
fact, they are sensitive to both mean and variance changes. The m-control distribution of 
i/ is N (0,1) 

The two standardized Cusum scale chart is defined as 

S* = max[o, v, - k + S* A ] 2 5 2 3.1 2 

S~ = max[o,-& ~v l +5'",] 2 523 1 3 

where Sq = S 0 ~ = 0 


2.6 MONITORING PROCESS VARIABILITY IN MULTIVARIATE 
PROCESSES : 

Two methods are also analyzed for monitoring process variability 

• Method Based on Sample Generalized Variance 

• Univariate S 2 Control Chart 
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2.6 1 Sample Generalized Variance Method: 


This method is used for monitoring multivariate process variability In multivariate 
piocess, process variability is summarized by the p x p covariance matrix 2, which is 
estimated by sample covariance matrix S, which is given by 


5 


2 

1 


s = 



S\ 






2611 


Sample covariance matrix S, is calculated by analysis of m-control preliminary data 
This method is based on the sample covariance matrix |S| 


Control Limits for Chart: 

UCL = |2|(A, +3b'( 1 ) 

Center line = b,|S| 2 6 12 

LCL = |2|(&,-3^) 


If we estimate 2 from sample covariance matrix S, then we take |S'|/Z> 1 as estimator of 
2 . b { and b 2 are given by 


4 > 


2 6.1 3 


(*-iy ,-i 
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1 


2.6 1 4 


b 2 = 


(n-l) 2p 


n(«-o 


7=1 7=1 J 


We calculate the sample covariance matrix for each sample similar to Eq-2.6 1 1 and 
values of is plotted for each sample on the control chart with above mentioned 
control limits 

This method shows drawbacks, when determinant of matrices are equal, even though 
elements of matrices are different i e they have different process variability and 
con elation between the variables Hence, it is wise to use univariate S 2 control charts m 
conjunction with this chart 
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CHAPTER 3 


PERFORMANCE EVALUATION 

This chapter is organized into two sections In the first section, computational performance of all charts 
are discussed which is based on the output obtained from the respective charts code. In the second 
section design of experiments are discussed for all chart Since the performance of charts seems to 
depend on the con elation coefficient, sample size and deviation in the process mean or process 
deviation, factorial design of experiment are performed for all charts and their performance at different 
levels of these three factors are analyzed 


3.1 Computational Performance 

Many problems m industrial quality control involve a vector of measurements of several characteristics 
simultaneously rather than a single characteristics For example, we may want to simultaneously control 
both the length and the diameter of the pipe In other words, both the length and diameter must be 
acceptable for the pipe to be usable Although one could monitor the process by separate control chart 
foi each quality characteristics, but results of monitoring depend on the correlation between 
chaiacten sties If the correlation is high, the result may not be acceptable. So, controlling each 
characteristic separately may not yield a product in which all variables should be acceptable In many 
industiies, they use separate control chart for each variable due to computational intricacies of 
multivariate control charts But, now with modern computers, one can easily use several types of 
multivariate charts The performance of different charts (some of which are multivariate charts and some 
are univariate charts) are computationally evaluated at different correlation coefficient, sample-size and 
deviations levels Umvauate charts are assumed not to give good results, when there is high correlation 
between variables So, they are analyzed at different correlation coefficient to observe the relationship. ; 
The performances of different multivariate charts are also evaluated to see the effect of three parameter 
(Con elation coefficient, sample size and deviations) 


Following charts for monitoring mean changes in the multivariate processes are analyzed: 


. CHI-SQUARE (z 2 ) CONTROL CHART 
. T_SQUARE (T 2 ) CONTROL CHART 
. MULTIVARIATE CUSUM CHARTS : 

(1) CUMULATIVE SUM OF T 

(2) CUSUM VECTOR SCHEME 

. MULTIVARIATE EXPONENTIALLY MOVING-AVERAGE CHART 
. REGRESSION-ADJUSTMENT PROCEDURE 

(1) Z-METHOD FOR MONITORING PROCESS MEAN 

(2) Y-METHOD 

. UNIVARIATE SHEWHART X-BAR 

. UNIVARIATE X-CUSUM CHART FOR MONITORING PROCESS MEAN 


Following charts for monitoring process variance changes in the multivariate processes are 
analyzed: 


• REGRESSION-ADJUSTMENT PROCEDURE 

1. Z-METHOD FOR MONITORING PROCESS VARIANCE 

• UNIVARIATE X RANGE CHART 

• MULTIVARIATE PROCESS VARIANCE CHART 

• UNIVARIATE 5 2 CONTROL CHART 

• UNIVARIATE X-CUSUM CHART FOR MONITORING PROCESS VARIANCE 


Codes are developed for all charts. Since original data (measurements of quality characteristics of 
multivariate process) are not available, so test data are generated from multivariate normal random 
generator The subroutine for multivariate normal random generator is available in the NAG library 
(written in FORTRAN language) Two sets of problems are designed, one for two variables and another 
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for five variables Perfoimance of these charts is evaluated for different values of correlation coefficient, 
sample size and deviation in process mean or variance 

Values of Factors : 

Following values foi different factors are selected 

Values of Correlation Coefficient Two Variables case: 

• Low value 

1.000000 0.100000 
0.100000 1.000000 

• Medium value 

1.000000 0.723009 
0.723009 1.000000 

• High value 

1.000000 0.900000 
0.900000 1.000000 


Values of Correlation Coefficient for Five Variables case : 

• Low value 

1.000000 0.100000 0.100000 0.100000 0.100000 

0.100000 1.000000 0.100000 0.100000 0.100000 

0.100000 0.100000 1.000000 0.100000 0.100000 

0.100000 0.100000 0.100000 1.000000 0.100000 

0.100000 0.100000 0.100000 0.100000 1.000000 

• Mixed value 

1.000000 0.138800 0.349600 0.082900 0.265200 

0.138800 1.000000 0.732400 0.913000 0.693200 

0.349600 0.732400 1.000000 0.682400 0.821400 

0.082900 0.913000 0.682400 1.000000 0.764000 

0.265200 0.693200 0.821400 0.764000 1.000000 


• High value 


1.000000 

0.900000 

0.900000 

0.900000 

0.900000 

0.900000 

1.000000 

0.900000 

0.900000 

0.900000 

0.900000 

0.900000 

1.000000 

0.900000 

0.900000 

0.900000 

0.900000 

0.900000 

1.000000 

0.900000 

0.900000 

0.900000 

0.900000 

0.900000 

1.000000 


Sample Sizes : 

Low (1), medium (6) and High (10) 

Deviations : 

Deviation for process means and variance are measured as multiple of process 
variance. Combinations of 0.25, 0.5, 1.0, 1.5, and 2.0 of process variance (<x) ) both positive and 
negative were used to generate deviation from the mean and variance. 


Initially, using NAG library subroutine sixty base data are generated, which are m-control for all charts 
Then, specified deviation is given after 50 th sample and signaling capabilities of different charts are 
analyzed The charts are peimitted to rim for 10 more samples after giving deviations Standard values 
of means and covariance matrix are used The charts for monitoring mean deviation are also analyzed 
when the variation is made in the process variance Similarly charts for monitoring process variance 
deviation are also analyzed when the deviation is made in the process means 


Table- 1 gives total number of problems m which process mean deviations were introduced for two ana 
five number of variables 


Table 3.1 NUMBER OF PROBLEMS FOR MEAN CHANGES 


Number of Variables 

2 

5 

Total No. of Problems 

46 * 9 

64 *9 

Total No. of Positive Problems 

17 *9 

27 * 9 

Total No. of Negative Problems 

17*9 

27 * 9 

Total No. of Mixed Problems 

12 *9 

10*9 


Table-3 2 to Table-3 5 represents computational output (Total number of signals out of total number o 
problems) for various charts for different types of problems 
There are three types of problems for mean deviations 

• In which only positive deviations were introduced m all variables 

• In which only negative deviation were introduced in all variables 

• In which both positive and negative deviations were introduced 


Table 3.2 OUTPUT (MEAN CHANGES) FOR NUMBER OF VARIABLE-2 


Cor- 

Sample 

T2_Chart 

^ 2 _Chart 

Cusum #1 

Cusum #2 

coff. 

Size(n) 

+ 

- 

± 

+ 

- 

± 

+ 

- 

± 

+ 

- 

± 

0 10 

1 

14 

0 

7 

12 

3 

7 

17 

10 

10 

15 

7 

9 


6 

16 

16 

12 

16 

15 

12 

17 

15 

12 

17 

15 

12 


10 

16 

16 

12 

17 

16 

12 

17 

16 

12 

17 

16 

12 

0 723 

1 

17 

8 

11 

16 

6 

10 

17 

11 

12 

16 

10 

12 


6 

17 

17 

12 

17 

17 

12 

17 

16 

12 

17 

17 

12 


10 

17 

17 

12 

17 

17 

12 

17 

17 

12 

17 

17 

12 

0 90 

1 

17 

15 

12 

16 

11 

12 

17 

13 

12 

16 

12 

12 


6 

17 

17 

12 

17 

17 

12 

17 

17 

12 

17 

17 

12 


10 

17 

17 

12 

17 

17 

12 

17 

17 

12 

17 

17 

12 
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Table 3.3 OUTPUT (MEAN CHANGES) FOR NUMBER OF VARIABLE-2 


Corr- 

Coff. 

Smaple 

Size(n) 

Expon_Chart 

ZMethod 

XJBar Chart 

X_Cusum Chari 

+ 

- 

± 

+ 

- 

± 

+ 

- 

± 

+ 

- 

± 

0 10 

1 

17 

10 

10 

16 

10 

10 

15 

9 

9 

17 

10 

10 


6 

15 

7 

10 

17 

15 

12 

17 

15 

12 

17 

15 

12 


10 

14 

12 

10 

17 

16 

12 

_ 

17 

16 

12 

17 

16 

12 

0 723 

1 

17 

8 

12 

16 

7 

12 

16 

6 

11 

16 

6 

10 


6 

15 

12 

12 

17 

17 

12 

17 

16 

12 

17 

16 

12 


10 

14 

12 

12 

16 

17 

12 

16 

17 

12 

17 

16 

12 

0 90 

1 

17 

10 

12 

13 

13 

12 

17 

17 

12 

17 

7 

10 


6 

15 

14 

12 

17 

17 

12 

17 

13 

12 

17 

13 

12 


10 

15 

15 

12 

17 

17 

12 

17 

16 

12 

17 

16 

12 


Table 3.4 OUTPUT (MEAN CHANGES) FOR NUMBER OF VARIABLE-5 


Cor- 

Sample 

T2_Chart 

X 2 JC hart 

Cusum #1 

Cusum #2 

Coff. 

Size(n) 

+ 

- 

± 

+ 

- 

± 

+ 

- 

+ 

+ 

- 

± 

0 10 

1 

12 

6 

9 

17 

6 

9 

24 

17 

9 

26 

18 

9 


6 

23 

22 

10 

24 

23 

10 

26 

26 

10 

24 

24 

10 


10 

23 

24" 

10 

24 

24 

10 

26 

27 

10 

26 

27 

10 

0 90 

1 

24 

23 

10 

27 

23 

10 

27 

23 

10 

27 

27 

10 


6 

23 

24 

i 

10 

25 

25 

10 

24 

27 

10 

27 

27 

10 


10 

24 

26 

10 

24 

27 

10 

24 

27 

10 

27 

27 

10 
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Table 3.5 OUTPUT (MEAN CHANGES) FOR NUMBER OF VARIABLE-5 


Corr 

Coff. 

Sample 

Size(n) 

Z 

+ 

_Meth< 

id 

± 

x_ 

+ 

Bar Ch 

lart 

± 

x_c 

+ 

usum C 

lhart 

± 

0 10 

1 

25 

19 

9 

17 

19 

8 

25 

19 

9 

6 

26 

24 

10 

26 

22 

10 

26 

27 

10 

10 

26 

27 

10 

25 

26 

10 

26 

27 

10 

0 90 

1 

20 

10 

10 

26 

18 

10 

26 

17 

10 

6 

20 

18 

10 

24 

23 

10 

26 

27 

10 

10 

20 

20 

10 

26 

24 

10 

26 

24 

10 


Z_Method for momtoung process mean deviation, univariate Shewhart X_Bar chart and umvariati 
X_Cusun chart also provide the idea about out-of-control and m-control variables However, it is no 
always true that these charts provide true result about out-of-control and m-control variables The result: 
seem to depend on correlation coefficient, sample size and magnitude of deviation introduced m th< 
process mean Table-6 provides the number of problems in which these charts gave true result fo 
different magnitude of correlation coefficient, sample size and number of variable 
Table 3.6 Output (Mean changes) for True Differentiation Between 


Out-of-Control and In-Control Variables (Numbar of Variable = 2) 


Cor- 

Sample 

Zjmethod 

X_Bar Chart 

X_Cusum 

Coff. 

Size(n) 

+ 

- 

± 

+ 

- 

± 

+ 

- 

± 

0 10 


12 

4 

9 

9 

3 

7 

13 

5 

8 


6 

14 

13 

12 

17 

14 

11 

17 

14 

11 


10 

14 

15 

12 

14 

15 

11 

17 

16 

12 

0 723 

1 

12 

12 

12 

17 

3 

5 

16 

3 

5 


6 

12 

12 

12 

16 

15 

11 

16 

15 

11 


10 

12 

12 

12 

17 

15 

12 

17 

16 

12 

0 90 

1 

12 

12 

12 

1 

17 

15 

12 

17 

4 

5 


6 

12 

12 

12 

17 

11 

10 

17 

11 

10 


10 

12 

12 

12 

17 

15 

11 

17 

16 

11 




Table 3.7 Output (Mean changes) for True Differentiation Between 


Out-of-Control and In-Control Variables (Numbar of Variable = 5) 


Cor- 

Sample 

z 

T3 

O 

<u 

B 

i 

X_Bar Chart 

X Cusum Chart 

Coff. 

Size 

+ 

- 

± 

+ 

- 

+ 

+ 

- 

± 

0 10 

1 

16 

4 

0 

4 

1 

0 

20 

10 

2 


6 

15 

12 

7 

20 

13 

3 

22 

25 

8 


10 

15 

15 

8 

13 

14 

4 

25 

22 

10 

0 90 

1 

12 

12 

12 

11 

8 

0 

25 

8 

5 


6 

12 

12 

12 

16 

13 

4 

22 

22 

10 


10 

12 

12 

12 

18 

25 

8 

24 

27 

6 


Z_Method for monitoring process variance deviation, univariate X_Range method, univariate 
X_Cusum and univariate S 2 control chart are used for monitoring deviation introduced in the process 
variance These charts suppose to give true differentiation between out-of-control and m-control 
variables But for these charts also as in charts for monitoring process mean deviation, it does not hold 
all the time that these provide true results 

Table-3 8 gives total number of problems in which process mean deviations were introduced for two and 
five number of variables 

Table 3.8 NUMBER OF PROBLEMS FOR VARIANCE CHANGES 


Number of Variables 

2 

5 

Total No. of Problems 

22* 9 

44 *9 


Table-3 9 and Table-3 10 represents computational output (Total number of signals out of total number 
of problems) for various charts for different types of problems Only positive deviations are introduced 
m the process variance 
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Table 3.9 


OUTPUT FOR NUMBER OF VARIABLE - 2 



Table 3.10 OUTPUT FOR NUMBER OF VARIABLE - 5 




10 

6 

34 

14 

32 

0 90 

1 

6 

32 

27 



6 

6 

33 

31 

33 


10 

6 

34 

26 

34 





3.1 CHI-SQUARE ( Z 2 ) CONTROL CHART: 


• This chart gives almost satisfactory result at all conditions But, it gives better result at positrv 
deviations from target means than negative deviations especially at n = 1 as shown m Table 3.2 an 
Table-3 4 

• When theie is both positive and negative deviation in different variables, signals improve 

• For the same values of deviation, positive deviation is detected faster than negative deviation 

• When sample size and correlation coefficient increases detection in case of negative deviatio 
improves as shown in Table -3 2 and Table-3 4 

• It gives fake signal of process mean being out-of-control when deviation is introduced m the proces 
variance almost in all problems However there is difference m the magnitude Z 2 statistics. In cas 
of mean deviation, the magnitude of statistics increases continuously after deviation if it is run afte 
introducing deviation But, m the case of variance deviation, its magnitude increases just afte 
deviation After that all values do not fall out-of-control if it is run after deviation Only some of th 
values fall out-of-conti ol and others remain m control So there is no continuous increment o 
magnitude On this basis deviation type may be recognized 

3.2 T SQUARE (T 2 ) CONTROL CHART : 

• T 2 chart is similar to Chi-Square chart, except that estimated value of means and covariance matn: 
from past in-control data rather than standard value of these are used 

• Results are almost same as in Chi-Square chart It also gives better result at positive deviations fron 
target means than negative deviations especially at n = 1 as shown m Table 3 2 and Table-3 4 

• When there is both positive and negative deviation m variables, improved signal is observed. 

• For the same values of deviation, positive deviations are detected faster than negative deviation 

• When sample size and correlation coefficient increase detection of negative deviation also improves 
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As in Chi_Square control chart, it also gives fake signal of process mean being out-of-control when 
deviation is introduced in the process variance almost m all problems However, in case of mean 
deviation, the magnitude of T~ statistics increases continuously if it is run after introducing mean 
deviation But, m the case of variance deviation, its magnitude increases just after deviation After 
that all values do not fall out-of-conti ol if it is run after deviation Only some of the values fall out- 
of-control and others remain in control So there is no continuous increment of magnitude On this 
basis deviation type may be lecognized 


3.3 MULTIVARIATE CUSUM CHARTS : 


3.3.1 CUMULATIVE SUM OF T AND CUSUM VECTOR SCHEME : 


• Cusum Vector method gives faster signals than Cumulative Sum of T (COT) atn=l 

• As sample size n increases, COT scheme starts to give faster signal than Cusum Vector method. 

• It gives fake signal of process mean being out-of-control when deviation is introduced in the process 
variance 

• At n = 1, Cusum vector method did not give signal of process being out-of- control even after 10 
runs aftei process deviation was introduced in more problems than the problems in which it gives 
such signal 

• As n increases, both starts to give signal, but COT signals are faster 

• At n = 1, these charts performs better than other chaits However as correlation coefficient increases 
they perform similar to other charts as shown in Table 3 2 and Table-3 4 

• It gives fake signal of process mean being out-of-control when deviation is introduced in the process 
vaiiance 

• These chaits signal in more test problems with negative deviation than T 2 chart at n = 1. As 
correlation coefficient increases T2 chart signals better. 


49 



3.4 MULTIVARIATE EXPONENTIALLY MOVING-AVERAGE CHART : 


• This chart pei forms comparable to others charts at n — 1 as shown m Table-3 3, but as n increases it 
signal are not as fast as of other charts 

• With the increase of con-elation coefficient, the signals are faster 

• This chart signals in moie test problems with negative deviation than T 2 chart at n = 1 As 
conelation coefficient increases T2 chart signals better 

• It gives fake signal of process mean being out-of-control when deviation is introduced m the process 
variance 

3.4 REGRESSION-ADJUSTMENT PROCEDURE : 

3.5.1 Z-METHOD FOR MONITORING MEAN DEVIATION: 

• This method is good m the sense that, separate diagnostic procedure is needed this method This 
chait gives clear identification of out-of-control variable or set of variables. 

• But, this does not hold all the time When conelation is low i e less than 0 5, we get good 
differentiation between out-of-control and in-control variables As conelation increases, out-of- 
control variables also affect other variables 

• At low correlation coefficient this chart performs similar to Shewhart X_Bar chart, but as conelation 
coefficient increases, deviation diffused in other variables also which does not happen in case of 
X_Bai chart and univariate X_Cusum chart as shown in Table-3 6 and Table-3.7 At large number of 
variable also, this diffusion occurs in this chart Hence at high conelation coefficient and large 
numbei of variables this chart does not perform better as X_Bar and X_Cusum performs 

• It gives fake signal of process mean being out-of-control when deviation is introduced m the process 
variance almost in all cases Since statistics of this chart is based on cumulative sum Hence there is 
continuous increment of statistics even in the case of process variance deviation as in the case of 
process mean deviation. 
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3.5.2 UNIVARIATE SHEWHART X-BAR CHART : 


• Separate univariate X_bar chart is run for all p variables to observe the difference in the result due 
to correlation coefficient 

• This chait gives more cleai result about which variable or set of variables are out-of-control 

• It gives fake signal of process mean being out-of-control when deviation is introduced in the process 
variance almost m all cases 


3.5.3 UNIVARIATE X-CUSUM CHART : 


• Separate Univariate X_Cusum charts are also maintained for all p variables 

• This chait performs satisfactory even at 0 9 correlation coefficient as shown in Table-3 6 and 
Table-3 7 

• This chart better than Z Method and X Bar chart as it is clear from Table-3 6 and Table-3 7. 


In the next section, we discuss the perfoimance of charts designed to detect process variance deviation 


3.6 MONITORING MULTIVARIATE PROCESS VARIANCE 


3.6.0 Z-METHOD FOR MONITORING PROCESS VARIANCE 


• At low coirelation coefficient and number of variables, there is more number of clear differentiation 

between out-of-control and in-control variables, but as they increases, deviation diffused into other 

variables also as shown in Table-3 9 and Table-3 10. So there is less number of clear signals Other 

charts perform better than this chart. tufcN f KAL. LiBHAB) 

S. I. T», KAMWW 

•m, 
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3.6.1 MULTIVARIATE PROCESS VARIANCE CHART : 


• It gives fake signal of process variance being out-of-control when deviation is introduced in the 
process mean almost in all cases 

• This chart also does not identify which variable or set of variables is out-of-control It gives only 
oveiall mfoimation that process is going out-of-control 

• When this chart is mn for process variance deviation, it gives very good result After one run after 
deviation it gives out-of-conti ol signal 


3.6.2 UNIVARIATE S 2 CONTROL CHART : 


• The drawback of above mentioned variance chart is described in Section-2 6 1 It’s suggested to use 
above mentioned variance chart in conjunction of univariate S 2 control chart 

• This chart gives better result than sample covariance method 

• It has diagnostic capability, to identify variable or set of variables which ls/are out-of-control. 

• This chart perform satisfactoiy all correlation and sample size but it perform better al low number of 
variables as there is large number of clear signal of out-of-control variables as shown m Table-3.9 
and Table-3 10 

• When this chart is run for process mean deviation, out-of-control signals are generally not observed 

• This chart performs better than Z_Method and univariate X_Cusum chart 
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3.6.3 UNIVARIATE X-RANGE CHART : 


• Separate univariate X-Range chait is mn for each of p variables 

• This chart does not perform as good as S' 2 control chait when process mean deviation is introduced 
This gives out-of-control signal when deviation is introduced in process mean 

• This chart performs similar to S 2 control chart in variance deviation as shown m Table-3 9 and 
Table-3 10 
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3.7 Design of Experiment 


As discussed in Section-3.1, all charts are permitted to run for different combinations of 
factors The results are tabulated With these observations a factorial design of 
experiment (DOE) with three factors correlation coefficient, sample-size and deviations 
was conducted The details of factor levels as discussed above, also shown in Table-3 1 


Table 3.11 Factors of Design of Experiment 


Factors 

Levels 

Correlation Coefficient 

Low, Medium and High 

Sample size 

1, 6 and 10 

Deviation 

20 Combinations of Low, Medium and 
High 


Initially DOE is conducted for two variable cases The output are taken m the form of +1 
(if a chart signals after introducing deviation) and -1 (if a chart does not signal after 
intioducing deviation) Each expenment is run for 10 additional samples after the 
deviation is introduced DOE is conducted only for main effects Thus total of 180 
(3 x 3 x 20) problems were solved for different charts Table-3 2 to Table-3 10 show the 
summary of effects of the experimental design for different charts 


3.8 Analysis of DOE 

From the design of experiment, following results are obtained for different charts 

3.8.1 T 2 CONTROL CHART: 

From the design of experiment it can be observed that following factors are significant 

• Correlation Coefficient 

• Sample Size 
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Table 3.12 Summary of Effects for T 2 Control Chart 


Effect 

SS 

Degr of 
Freedom 

MS 

F 

P 

Intercept 

108 8889 

1 

108 8889 

391 1975 

0 

COR COF 

10 97778 

2 

5 488889 

19 71955 

2 32E-08 

SAM SIZE 

8711111 

2 

4 355556 

15 6479 

6 41E-07 

DEVI ATI 0 

8 

19 

0 421053 

1 512687 

0 08773 

Error 

43 42222 

156 

0 278348 




Following aie the obseivations made from the results as shown in Fig-3 1 to Fig-3 3 

With the increase in correlation coefficient, this chart performs better i.e. signal 
for out-of-control becomes faster. 

As sample size n increases, performance increases. But after sample size n = 6, 
improvement m performance is constant. 

At high deviation, signals are faster than low deviation. 


COR_COF, LS Means 
Current effect F(2, 156)=19 720, p= 00000 
Effective hypothesis decomposition 


ce 

< 

x 

o 


1 3 
1 2 
1 1 
1 0 
09 
0 8 
07 
0 6 
05 
04 
03 

02 



COR_COF 


Figure 3.1 T 2 CHART vs. CORRELATION COEFFICIENT 
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T2_CHART ^ T2_CHART 


SAM_S1ZE, LS Means 
Current effect F(2, 156)=15 648, p= 00000 
Effective hypothesis decomposition 



SAM SIZE 


igure 3.2 


T 2 CHART vs. SAMPLE SIZE 


DEVIATIO, LS Means 

Current effect F(19, 156)=1 5127, p= 08773 
Effective hypothesis decomposition 



4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 


Figure 3.3 


T 2 CHART vs. DEVIATION 





3.8.2 MULTIVARIATE CUSUM CHARTS 


3.8.1.1 CUSUM #1 : 


Fiom the design of experiment it can be observed that following factors are significant. 


• Deviation 

• Sample Size 


Table 3.13 Summary of Effects for Cusurn#! 


Effect 

SS 

Degr of 
Freedom 

MS 

F 

P 

Intercept 

131 7556 

1 

131 7556 

672 1831 

0 

COR COF 

0 844444 

2 

0 422222 

2 15407 

0 119449 

SAM SIZE 

3 244444 

2 

1.622222 

8 276163 

0 000384 

DEVI ATI 0 

13 57778 

19 

0 71462 

3 645808 

3 45E-06 

Error 

30 57778 

156 

0 196011 




Following are the observations made fiom the lesults as shown in Fig-3 4 to Fig-3 6 

• With the increase in correlation coefficient, this chart performs better. 

• As sample size n increases, performance increases. After n = 6, performance 
improves slowly. 

• At high deviation, signals are faster than low deviation. 
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yvHo ion v* who ion 


COR_COF, LS Means 
Current effect F(2, 156)=2 1541, p= 11945 
Effective hypothesis decomposition 



gure 3.4 CUSUM#1 vs. CORRELATION COEFFICIENT 


SAM_S1ZE, LS Means 
Current effect F(2, 156)=8 2762, p= 00038 
Effective hypothesis decomposition 



SAM SIZE 


Figure 3.5 


CUSUM#I vs. SAMPLE SIZE 





DEVIATIO, LS Means 

Current effect F(19, 156)=3 6458, p= 00000 
Effective hypothesis decomposition 



Figure 3.6 CUSUM#1 vs. DEVIATION 


S.8.2.2 CUSUM #2: 

From the design of experiment it can be observed that following factors are significant 

• Correlation 

• Deviation 

• Sample Size 


Table 3.14 Summary of Effects for Cusum #2 


Effect 

SS 

Degr of 
Freedom 

MS 

F 

_2 

Intercept 

121 6889 

1 

121 6889 

513 9928 

0 

COR COF 

4 444444 

2 

2 222222 

9 386282 

0 000142 

SAM SIZE 

7 511111 

2 

3 755556 

15 86282 

5 36E-07 

DEVIATION 

9 422222 

19 

0 495906 

2 094623 

0 007231 

Error 

36 93333 

156 

0 236752 




Following are the observations made from the results as shown m Fig-3 7 to Fig-3 9 

• With the increase in correlation coefficient, this chart performs better. 

• After 0.723 correlation coefficient , improvement in performance is constant. 

• As sample size n increases, performance increases. After n = 6, improvement is 
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MC2_CHAR go MC 2 _CHAR 


constant. 

At high deviation, signals are faster than low deviation. 


COR.COF, LS Means 
Current effect F(2, 1 56)=9 3863, p= 00014 
Effective hypothesis decomposition 



;ure 3.7 


CUSUM#2 vs. CORRELATION COEFFICIENT 


SAM^SIZE, LS Means 
Current effect F(2, 156)=15 863, p= 00000 
Effective hypothesis decomposition 



SAM SIZE 


CUSUM#2 vs. SAMPLE SIZE 


Figure 3.8 


m 





DEVIATIO, LS Means 

Current effect F(19, 156)=2 0946, p= 00723 
Effective hypothesis decomposition 



DEVIATIO 

Figure 3.9 CUSUM#2 vs. DEVIATION 

3.8.3 SHEWHART XBAR CHART : 

From the design of experiment it can be observed that following factors are 
significant 

• Deviation 

• Sample Size 


Table 3.15 Summary of Effects for X_Bar Chart 


Effect 

SS 

Degr of 
Freedom 

MS 

F 

P 

Intercept 

105 8 

1 

105 8 

437 4064 

0 

COR COF 

1 2 

2 

06 

2 480565 

0 086992 

SAM SIZE 

3 733333 

2 

1 866667 

7 717314 

0 000637 

DEVIATIO 

31 53333 

19 

1 659649 

6 861447 

6 31E-13 

Error 

37 73333 

156 

0 24188 




Following are the observations made from the results as shown m Fig-3 10 to Fig- 3 12 

• With the increase in correlation coefficient, this chart performs better. 

• As sample size n increases, improvement in performance is slower. 

• At high deviation, signals are faster than low deviation 
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DEVIATIO, LS Means 

Current effect F{19, 156)=6 8614, p= 00000 
Effective hypothesis decomposition 



DEVIATIO 


Figure 3.12 XJBAR CHART vs. DEVIATION 


3.8.4 UNIVARIATE X_CUSUM CHART : 

From the design of expenment it can be observed that following factors are significant. 

• Deviation 

• Sample Size 


Table 3.16 Summary of Effects for X_Cusum Chart 


Effect 

SS 

Degr of 
Freedom 

MS 

F 

9 

Intercept ^ 

93 88889 

1 

93 88889 

433 0486 

0 

COR COF 

0 311111 

2 

0 15555CT 

0 717477 

0.489585 

SAM SIZE 

7 644444 

2 

3 822222 

17 62943 

1 25E-07 

DEVIATION 

44 33333 

19 1 

2 333333 

10 76216 

0 

Error 

33 82222 

156 

0 216809 




Following are the observations made from the results as shown in Fig-3 1 j to Fig-3.15 

• With the increase in correlation coefficient, this chart performs better. 

• As sample size n increases, improvement in performance is slower. 

• At high deviation, signals are faster than low deviation 
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wnsno x 


COR_COF, LS Means 
Current effect F(2, 156)= 71748, p= 48958 
Effective hypothesis decomposition 



Figure 3.13 X^CUSUM vs. CORRELATION COEFFICIENT 


SAM_SIZE, LS Means 
Current effect F(2, 156)=17 629, p= 00000 
Effective hypothesis decomposition 



1 2 3 

SAM SIZE 


Figure 3.14 


X CUSUM vs. SAMPLE SIZE 





DEVIATiO, LS Means 

Current effect F(19, 156)=10 762, p=0 0000 
Effective hypothesis decomposition 



DEVIATiO 


Figure 3.15 X_CUSUM vs. DEVIATION 


3.8.5 ZMETHOD : 


From the design of experiment it can be observed that following factors are significant- 

• Deviation 

• Sample Size 


Table 3.17 Summary of Effects for Z_Method 


Effect 

SS 

Degr of 
Freedom 

MS 

F 

P 

Intercept 

125 

1 

125 

543.0074 

0 

COR COF 

0 933333 

2 

0 466667 

2.027228 

0 135155 

SAM SIZE 

4 933333 

2 

2 466667 

10 71535 

4 36E-05 

DEVIATION 

13 22222 

19 

0 695906 

3 023059 

7 94E-05 

Error 

35 91111 

156 

0 230199 




Following are the observations made from the results as shown in Fig-3 16 to Fig-3.18 

• With the increase in correlation coefficient, this chart performs better. 

• As sample size n increases, improvement in performance starts to decrease. 
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Z_METHOD 


At high deviation, signals are faster than low deviation 


1 1 


i o 


0 9 

a 

o 

x o 8 

LU 



06 

05 


COR_COF, LS Means 
Current effect F(2, 156)=2 0272, p= 13516 
Effective hypothesis decomposition 



1 2 3 

CORJSOF 


Figure 3.16 ZJMETHOD vs. CORRELATION COEFFICIENT 


SAM^SJZE, LS Means 
Current effect F(2, 156)=10 715, p= 00004 
Effective hypothesis decomposition 


1 2 
1 1 
1 0 
0 9 
0 8 
I 07 
0 6 
05 
04 



1 


2 

SAM_SIZE 


3 


Figure 3.17 Z METHOD vs. SAMPLE SIZE 





DE VIATIO, LS Means 

Current effect F(19, 156)=3 0231, p= 00008 
Effective hypothesis decomposition 



DE VIATIO 

Figure 3.18 ZMETHOD vs. DEVIATION 

3.8.6 EXP ON CH ART : 

From the design of experiment it can be observed that following factors are significant 

• Deviation 


Table 3.18 Summary of Effects for ExponChart 


Effect 

SS 

Degr of 
Freedom 

MS 

F 

P 

Intercept 

1 1 8642 

1 

1 1 8642 

29 81204 

7 26E-07 

COR COF 

1 876543 

2 

0 938272 

2 357664 

0 102334 

SAM SIZE 

0 395062 

2 

0 197531 

0 49635 

0 610937 

DEVIATIO 

39 80247 

8 

4 975309 

12 50182 

7 65E-11 

Error 

27 06173 

68 

0 397967 




Following are the observations made from the results as shown m Figure 3 19 to 3 21 

• With the increase in correlation coefficient, this chart performs better. 

• As sample size n increases, improvement starts to decrease and after n = 6 it 
becomes constant. 

• At high deviation, signals are faster than low deviation 
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3.9 Design of Experiment with Detection ( Nondetection ) as Performance Measure 


From the previous DOE, it is observed that deviation m fact can be clustered into 3 or 4 
levels. So in this set of experiment, we select the detection of signal as two-performance 
level For this purpose if signal is detection 10 samples after the deviation is introduced, 
an output value of +1 (Detection) is assigned otherwise output value of -1 
(Nondetection) is assigned Also, this DOE considered all effects i e interactions effects 
also The DOE is replicated 5 times and analyzed for six charts Thus total of 180 
problems (3 * 3 * 4 * 5) were solved for each chart. 


Table 3.19 


Factors of Design of Experiment 


Factors 

Levels 

Correlation Coefficient 

Low, Medium and High 

Sample Size 

1, 6 and 10 

Deviation 

4 Combinations of Low, Medium and 
High 


Levels of Deviation: 

• Low < ± 0.5 

• Medium > ± 0.5 and < ± 1.0 

• High > ± 1.0 

• High and Low Combinations of Deviation 


3.10 Analysis of DOE 

From the design of experiment, following results are obtained for different charts 
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3.10.1 T 2 CONTROL CHART : 


From the design of experiment it can be observed that following factors are significant 

• Correlation Coefficient 

• Sample Size 

• Deviation 

• Correlation Coefficient and Sample Size 

• Correlation Coefficient and Deviation 


Table 3.20 Summary of Effects for T 2 Control Chart 


Effect 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 

3 622222 

144 

0.255556 

14 17391 

2 4E-06 

2 

2 

4 355556 

144 

0 255556 

17 04348 

2 27E-07 

3 

3 

1 481481 

144 

0 255556 

5 797101 

0 000906 

12 

4 

1 088889 

144 

0 255556 

4 26087 

0 00274 

13 

6 

1 103704 

144 

0 255556 

4 318841 

0 000489 

23 

6 

0 148148 

144 

0 255556 

0 57971 

0 746052 

123 

12 

0 17037 

144 

0 255556 

0 666667 

0 781 


Following are the observations made from the results as shown in Fig-3 22 to Fig- 3.26 

• With the increase in correlation coefficient, this chart performs better. 

• As sample size n increases, performance increases. But after sample size n = 6, 
improvement in performance is constant. 

• In the figure 3.25, third level of deviation factor is combination of high positive 
and high negative. T2 chart does not signal properly at negative high deviation, 
that’s why output is less than fourth level at which is T2 chart performs good. 

. In two-way interaction effect between sample size and correlation coefficient, at 
particular sample size, as correlation coefficient increases, performance 
improves. At n = 1, number of signals are more due to increase in correlation 
coefficient as compared to n > 1. 

. In two-way interaction effect between deviation and correlation coefficient, at 
particular deviation, as correlation coefficient increases, performance improves. 
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Variable T2_CHART Jg- Variable T2_CHART 


Plot of Means (unweighted) 
COR_COF Main Effect 
F(2,144)=14 17, p< 0000 



;ure 3.22 T 2 CHART vs. CORRELATION COEFFICIENT 


Plot of Means (unweighted) 
SAM_SIZE Mam Effect 
F(2,144)=17 04, p< 0000 



Figure 3.23 


T 2 CHART vs. SAMPLE SIZE 





Variable T2 CHART 


Plot of Means (unweighted) 
DEVIATIO Main Effect 
F(3,144)=5 80, p< 0009 



Figure 3.24 T 2 CHART vs. DEVIATION 


Plot of Means (unweighted) 
2-way interaction 
F(4,144)=4 26, p< 0027 



-o- CORJOOF 
G__l 1 

□ COR_COF 
G_2 2 

o CORJOQF 
G_3 3 


Figure 3.25 T 2 CHART vs. CORRELATION COEFFICIENT 

AND SAMPLE SIZE 





Plot of Means (unweighted) 
2-way interaction 
F(6,144)=4 32, p< 0005 



-o- 

COR COF 


G_1 1 

□ 

COR COF 


G_2 2 


COR COF 
G_3 3 


Figure 3.26 T 2 CHART vs. CORRELATION COEFFICIENT 

AND DEVIATION 


3.10.2 MULTIVARIATE CUSUM CHARTS 

3.10.2.1 CUSUM #1 : 

From the design of experiment it can be observed that following factors are significant 
» Sample Size 

• Deviation 

• Sample Size and Deviation 


Table 3.21 Summary of Effects for Cusum #1 


Effect 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 

0 288889 

144 

0.233333 

1 238095 

0 293003 

2 

2 

1 622222 

144 

0 233333 

6 952381 

0 001311 

3 

3 

1 622222 

144 

0 233333 

6 952381 

0 000211 

12 

4 

0 055556 

144 

0 233333 

0 238095 

0 916421 

13 

6 

0 288889 

144 

0 233333 

1 238095 

0 290286 

23 

6 

0 555556 

144 

0 233333 

2 380952 

0 031861 

123 

12 

0 055556 

144 

0 233333 

0 238095 

0 996041 
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Following are the observations made from the results as shown in Fig-3 27 to Fig-3 30 


. With the increase in correlation coefficient, this chart performs better. 

. As sample size n increases, performance increases. But after sample size n = 6, 
performance improvement is slow. 

. With increase in the deviation, charts signals are fast. 

• In the 2-way interaction effect between sample size and deviation, at low 
deviation as sample size increases signals are becoming faster, but as deviation 
increases increment of sample size does not much effect. 


Plot of Means (unweighted) 
CORJDOF Main Effect 
F(2,144)=1 24, p< 2930 



Figure 3.27 CUSUM #1 vs. CORRELATION COEFFICIENT 
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Variable MCI CHAR 


Plot of Means (unweighted) 
SAM_SIZE Mam Effect 
F(2,144)=6 95, p< 0013 



Figure 3.28 CUSUM #1 vs. SAMPLE SIZE 


Plot of Means (unweighted) 
DEVIATIO Mam Effect 
F(3,144)=6 95, p< 0002 



Figure 3.29 


CUSUM #1 vs DEVIATION 



Plot of Means (unweighted) 
2-way interaction 
F(6,144)=2 38, p< 0319 



-o- SAM_SIZE 
G_1 1 

o SAM_SIZE 
G_2 2 

o SAM__SIZE 
G_3 3 


Figure 3.30 CUSUM #1 vs SAMPLE SIZE AND DEVIATION 


3.10.2.2 CUSUM #2 

From the design of experiment it can be observed that following factors are significant: 

• Correlation Coefficient 

• Sample Size 

• Deviation 

• Sample Size and Deviation 

• Correlation Coefficient and Deviation 


Table 3.22 Summary of Effects for Cusum#2 


Effect 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 

08 

144 

02 

4 

0.020387 

2 

2 

3 2 

144 

02 

16 

5 31E-07 

3 

3 

2 214815 

144 

02 

1 1 07407 

1 38E-06 


4 

02 


02 

1 

0 409714 


6 

0 503704 


02 

2 518518 

0 023908 


6 

1 125926 


0.2 

5 62963 

2 82E-05 

123 

12 

0 081481 

144 

02 

0 407407 

0.958878 
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Following are the observations made from the results as shown m Fig-3 3 1 to Fig-3 35 

With the increase in correlation coefficient, this chart performs better. 

As sample size n increases, performance increases. But after sample size n = 6, 
performance improvement is slow. 

With increase in the deviation, charts’signals are fast. 

In two-way interaction effects between correlation coefficient and deviation, 
correlation coefficient has more effect than deviation. All correlation 
co efficient’ levels have the same effect. As both increases, they have same effect 
as they coincide. 

In two-way interaction effects between sample size and deviation, sample size 
has more effect than deviation. All sample size’levels have the same effect. As 
both increases, they have same effect as they coincide. 


Plot of Means (unweighted) 
COR_COF Mam Effect 
F(2,144)=4 00, p< 0204 



Figure 3.31 CUSUM#2 vs. CORRELATION COEFFICIENT 




Variable MC2 CHAR 


Plot of Means (unweighted) 
SAM_SIZE Main Effect 


F(2,144)=16 00, p< 0000 



Figure 3.32 CUSUM #2 vs. SAMPLE SIZE 


Plot of Means (unweighted) 
DEVIATIO Mam Effect 
F(3,144)=11 07, p<0000 



Figure 3.33 


CUSUM#2 vs. DEVIATION 
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3.10.3 Z METHOD : 


From the design of experiment it can be observed that following factors are significant 

• Sample Size 

• Deviation 

• Sample Size and Deviation 


Table 3.23 Summary of Effect for Z_Method 


Effect 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 

0 288889 

144 

0 233333 

1 238095 

0 293003 

2 

2 

2 955555 

144 

0 233333 

12 66667 

8 56E-06 

3 

3 

2 844445 

144 

0 233333 

12 19048 

3 73E-07 

12 

4 

0 055556 

144 

0 233333 

0 238095 

0 916421 

13 

6 

0 288889 

144 

0 233333 

1 238095 

0 290286 

23 

6 

1 177778 

144 

0 233333 

5 047619 

9 98E-05 

123 

12 

0 055556 

144 

0 233333 

0 238095 

0 996041 


Following are the observations made from the results as shown m Fig-3 36 to Fig-3 39 

• With the increase in correlation coefficient, this chart performs better. 

• As sample size n increases, performance deteriorate. 

• With increase in the deviation, charts ’signals are fast. 

• In two-way interaction effects between sample size and deviation, sample size 
has more effect than deviation. All sample size’levels have the same effect. As 
both increases, they have same effect as they coincide. 
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Variable Z_METHOD 


Plot of Means (unweighted] 
COR_COF Main Effect 


F(2,144)=1 24, p< 2930 



figure 3.36 


ZJVIETHOD 


vs. CORREL/* 


1 05 
1 00 
0 95 
0 90 
0 85 
0 80 
0 75 
0 70 
0 65 


Plot of Means (unweightec 
SAM_SIZE Main Effect 
F(2,144)=12 67, p< 0000 






Variable Z METHOD 


Plot of Means (unwei 
DEVI AT 10 Main Ef 
F(3,144)=12 19, p< 





3.10.4 SHEWHART XJ8AR CHART : 


From the design of experiment it can be observed that following factors are significant 


• Sample Size 

• Deviation 

• Correlation Coefficient and Sample Size 

• Correlation Coefficient, Sample Size and Deviation 


Table 3.24 Summary of Effects for X_Bar Chart 


Effect 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 

0 266667 

144 

0 233333 

1 142857 

0 321782 

2 

2 

0 866667 

144 

0 233333 

3 714286 

0 026738 

3 

3 

4 533333 

144 

0 233333 

19 42857 

1 24E-10 

12 

4 

0 933333 

144 

0 233333 

4 

0 004162 

13 

6 

0 266667 

144 

0 233333 

1 142857 

0 340595 

23 

6 

0 333333 

144 

0 233333 

1 428571 

0 207569 

123 

12 

0 666667 

144 

0 233333 

2 857143 

0 001488 


Following are the observations made from the results as shown in Fig-3.40 to Fig-3 44 

• With the increase in correlation coefficient, this chart performs better. 

• As sample size n increases, performance increases. 

• With increase in the deviation, charts ’signals are fast. 

• In two-way interaction effect between correlation coefficient and sample size, at 
low values both factor have same effect but as correlation coefficient increases, it 
has less effect than sample size. 

• In 3-way interaction effect among correlation coefficient, sample size and 
deviation, at high values all have same effect. 
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Variable X_BAR 


Plot of Means (unweighted) 
DEVIATIO Main Effect 
F(3,144)=19 43, p< 0000 



DEVIATIO 


Figure 3.42 


X BAR CHART vs. DEVIATION 


Plot of Means (unweighted) 
2-way interaction 
F(4 t 144)=4 00, p< 0042 



-O' COR COF 

GJ T 

a COR COF 
G_2 2 

o COR_COF 
G 33 


Figure 3.43 X BARCHART vs. CORRELATION COEFFICIENT 

AND SAMPLE SIZE 





Plot of Means 
3-way interaction 
F(12,144)=2 86, p< 0015 


£ 9 



G_1 1 G_3 3 G_1 1 G_3 3 G_1 1 G_3 3 G_1 1 G_3 3 

DEVIATIO DEVIATIO DEVIATIO DEVIATIO 

G_1 1 G_2 2 G_3 3 G_4 4 


-o- COR_COF 
G_1 1 

□ COR_COF 
G_2 2 

❖ COR_COF 
G_3 3 


Figure 3.44 X_BAR CHART vs. CORRELATION COEFFICIENT, 
DEVIATION AND SAMPLE SIZE 


3.10.5 UNIVARIATE X_CUSUM CHART : 

From the design of experiment it can be observed that following factors are significant 


• Sample Size 

• Deviation 

• Sample Size and Deviation 


Table 3.25 Summary of Effects for X_Cusum 


Effects 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 

0 266667 

144 

0 311111 

0 857143 

0 426526 

2 

2 

3 266667 

144 

0.311111 

105 

5 54E-05 

3 

3 

7 170371 

144 

0 311111 

23 04762 

3.04E-12 

12 

4 

0 133333 

144 

0 311111 

0 428571 

0 787819 

13 

6 

0 148148 

144 

0 311111 

0 47619 

0 825196 

23 

6 

1 192593 

144 

0 311111 

3.833333 

0 001413 

123 

12 

0.103704 

144 

0 311111 

0 333333 

0 981929 


87 







Following are the observations made from the results as shown m Fig-3 45 to Fig-3.48 


• With the increase in correlation coefficient, this chart performs worse. 

• As sample size n increases, performance also increases. 

• With increase in the deviation, charts’signals are fast. 

• In two-way interaction effects between sample size and deviation, sample size 
has more effect than deviation. As both factor increases, they have same effect as 
they coincide. 


Plot of Means (unweighted) 
COFLCOF Main Effect 
F(2 t 144)= 86, p< 4265 



Figure 3.45 X_CUSUM vs. CORRELATION COEFFICIENT 
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3.11 Design of Experiment with Run-Length as Performance Measure 

This design of experiment is done for two variable case, but run-length (number of 
samples after which there is signal) are considered as performance measure Levels of 
factors are same as m the DOE discussed in section-3.9 In case a signal in not observed 
m 10 samples, 11111 length is taken as 1 1 


Table 3.26 Factors of Design of Experiment 


Factors 

Levels 

Correlation Coefficient 

Low, Medium and High 

Sample Size 

1, 6 and 10 

Deviation 

4 Combinations of Low, Medium and 
High 


Levels of deviation: 

• Low < ± 0.5 

• Medium > ± 0.5 and < ± 1.0 

• High >±1.0 

• High and Low Combinations of Deviations 


3.12 Analysis of DOE 

Fiom the design of experiment, following results are obtained for different charts. 

3.12.1 T 2 CONTROL CHART : 

From the design of experiment it can be observed that following factors are significant 

• Correlation Coefficient 

• Sample Size 

• Deviation 

• Correlation Coefficient and Sample Size 
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Table 3.27 Summary of Effects for T 2 Control Chart 


Effect 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 

95 62222 

144 

4061111 

23 54583 

1 42E-09 

2 

2 

206 0722 

144 

4061111 

50 74282 

2 09E-17 


3 

118 1093 

144 

4 061111 

29 08299 

9 22E-15 


4 

22 38889 

144 

4 061111 

5.512996 

0 00037 

13 

6 

5 97037 

144 

4 061111 

1 470132 

0 192455 

23 

6 

2 42037 

144 

4 061111 

0 595987 

0.733166 

123 

12 

3 081481 

144 

4 061111 

0 758778 

0 691647 


Following are the observations made from the results as shown in Fig-3 49 to Fig-3 54 

• Graphs for this charts are similar to Figure from 3.22 to 3.24. which are 
conducted by using +1 and -1 as performance measure. 

• With the increase in correlation coefficient, this chart performs better. 

• As sample size n increases, performance increases. But after sample size n = 6, 
improvement in the performance is slow. 

• At high deviation, signals are faster than low deviation. 

. In two-way interaction effect between correlation coefficient and sample size, at 
particular sample size, as correlation coefficient increases, performance 
improves. At n = I, number of signals are more due to increase in correlation 
coefficient as compared to n > 1. 

• When DOE is conducted by using +1 and -I as performance measure, two-way 
interaction effect between correlation coefficient and deviation (section-3.10.1) 
was also significant, which is not here. 
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Variable T2 CHART 


Plot of Means (unweighted) 
DEVIATIO Mam Effect 
F(3,144)=29 08, p< 0000 


60 

5 5 

5 0 

fe 45 
< 

O 4 0 

P* 

^35 

0 ) 

JD 

S 30 
;> 

25 
20 
1 5 



G 1 1 


G 22 


G 33 


G 44 


DEVIATIO 


Figure 3.51 T 2 CHART vs. DEVIATION 


Plot of Means (unweighted) 
2-way interaction 
F(4,144)=5 51, p< 0004 



—o— COR COF, 

g_iT 

o CORJOOF 
G_2 2 

o COR COF 
G 33 


Figure 3.52 T 2 CHART vs. CORRELATION COEFFICIENT 

AND SAMPLE SIZE 





3.12.2 MULTIVARIATE CUSUM CONTROL CHARTS 


3.12.2.1 CUSUM #1 : 

Fiom the design of experiment it can be observed that following factors are significant. 

• Correlation Coefficient 

• Sample Size 

• Deviation 


Table 3.28 Summary of Effects for Cusum #1 Chart 


Effect 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 

23 43889 

144 

3 933333 

5 95904 

0 003263 

2 

2 

119 2389 

144 

3 933333 

30 31497 

1 03E-11 

3 

3 

135 1019 

144 

3 933333 

34 34793 

8.32E-17 

12 

4 

1 013889 

144 

3.933333 

0 257768 

0 904519 

13 

6 

1 97963 

144 

3.933333 

0 503296 

0 805098 

23 

6 

2 846296 

144 

3 933333 

0.723635 

0 631221 

123 

12 

0 465741 

144 

3 933333 

0 118409 

0 999886 


Following are the observations made from the results as shown in Fig-3 53 to Fig-3 55 

• With the increase in correlation coefficient, this chart performs better. 

• As sample size n increases, performance of this chart improves. But after sample 
size n = 6, improvement in performance is slow. 

• At high deviation, signals are faster than low deviation. 

• There are no significant interaction effects on this chart. 
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Variable MCI CHAR 


Plot of Means (unweighted) 
COR__COF Main Effect 
F(2,144)=5 96, p< 0033 



Figure 3.53 CUSUM#1 vs. CORRELATION COEFFICIENT 


Plot of Means (unweighted) 
SAM_SIZE Mam Effect 
F(2,144)=30 31 , p< 0000 



Figure 3.54 


CUSUM #1 vs. SAMPLE SIZE 




DEVIATIO Main Effect 



G_1 1 G_2 2 G_3 3 G_4 4 

DEVIATIO 

Figure 3.55 CL SUM #1 vs. DEVIATION 


3.12.2.2 CUSUM #2: 


Fiom the design of experiment it can be observed that following factors are significant 

• Correlation Coefficient 

• Sample Size 

• Deviation 


Table 3.29 Summary of Effects for Cusum#2 Chart 


Effect 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 

31 65 

144 

3 919445 

8 075124 

0 000474 

2 

2 

231 8 

144 

3 919445 

59 14103 

1 78E-19 

3 

3 

113 2889 

144 

3.919445 

28 90432 

1 09E-14 

12 

4 

7 925 

144 

3 919445 

2 02197 

0 094435 

13 

6 

2 205555 

144 

3 919445 

0 562721 

0 759412 

23 

6 

6 622222 

144 

3 919445 

1 689582 

0 127578 

123 

12 

1.347222 

144 

3 919445 

0 343728 

0.979414 
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Following are the observations made from the results as shown m Fig-3 56 to Fig-3 58 


• With the increase in correlation coefficient, this chart performs better. 

• As sample size n increases, performance increases. But after sample size n = 6, 
improvement in performance is slow. 

• At high deviation, signals are faster than low deviation. 


Plot of Means (unweighted) 
COFL.COF Mam Effect 
F(2,144)=8 08, p< 0005 



COR_COF 

Figure 3.56 CUSUM #2 vs. CORRELATION COEFFIECIENT 


Plot of Means (unweighted) 
SAM_SIZE Mam Effect 
F(2,144)=59 14, p< 0000 



Figure 3.57 Cusum #2 vs. Sample Size 
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DEVIATIO Mam Effect 
F(3,144)=28 90, p< 0000 



Figure 3.58 CUSUM#2 vs. DEVIATION 


3.12.3 ZMETHOD : 

Fiona the design of experiment it can be observed that following factors are significant 

• Correlation Coefficient 

• Sample Size 

• Deviation 


Table 3.30 Summary of Z_Method Chart 


Effect 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 

24 06667 

144 

3 755556 

6 408284 

0 002157 

2 

2 

134 3167 

144 

3 755556 

35 76479 

2 45E-13 

3 

3 

147.087 

144 

3 755556 

39.16519 

1 44E-18 

12 

4 

0 933333 

144 

3 755556 

0.248521 

0.910172 

13 

6 

3 503704 

144 

3 755556 

0.932939 

0 473365 

23 

6 

6 420371 

144 

3 755556 

1 709566 

0 12278 

123 

12 

1 57037 

144 

3 755556 

0.418146 

0 954494 


Following are the observations made from the results as shown in Fig-3 59 to Fig-3 61 


• With the increase in correlation coefficient, this chart performs better. 

• As sample size n increases, performance increases. But after sample size n = 6, 
improvement in performance is slow. 
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Plot of Means (unweighted) 
DEVIATIO Main Effect 
F(3,144)=39 17, p< 0000 



G_1 1 GJ2. 2 G_3 3 G_4 4 

DEVIATIO 


Figure 3.61 ZMETHOD vs. DEVIATION 

3.12.4 SHEWHART X BAR CHART : 


From the design of experiment it can be observed that following factors are significant. 

• Sample Size 

• Deviation 

• Correlation Coefficient and Sample Size 


Table 3.31 Summary of Effect of X_Bar Chart 


Effect 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 


2 

10 07222 

144 

3 755556 

2.681953 

0 071848 

HHH 

2 

112 1556 

144 

3 755556 

29 8639 

1 41E-11 

MEM 

3 

252 0352 

144 

3 755556 

67 10996 

3 29E-27 

12 

4 

24 80556 

144 

3.755556 

6 60503 

6.57E-05 

13 

6 

2 67963 

144 

3 755556 

0.713511 

0 639267 

23 

6 

6.540741 

144 

3 755556 

1 741617 

0 115428 

123 

12 

5 635185 

144 

3 755556 

1 500493 

0 13023 
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Following are the observations made from the results as shown in Fig-3 62 to Fig-3 65 

• With the increase in correlation coefficient, this chart performs better. 

• As sample size n increases, performance increases. But after sample size n = 6, 
improvement in performance is not much. 

• At high deviation, signals are faster than low deviation. 


Plot of Means (unweighted) 
COR__COF Main Effect 



G__1 1 G_2 2 G_3 3 

COR COF 


Figure 3.62 XJBAR CHART vs. CORRELATION COEFFICIENT 


Plot of Means (unweighted) 
SAM_SIZE Mam Effect 
F(2,144)=29 86, p< 0000 



Figure 3.63 XJBAR CHART vs. SAMPLE SIZE 
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Variable X BAR 


DEVI ATI O 


Figure 3.64 XJBAR CHART vs. DEVIATION 


Plot of Means (unweighted) 
2-way interaction 
F(4,144)=6 61, p< 0001 



~~o— CORJOOF 
G_1 1 

Q COR COR 
GJ2 2 

o COR_COF 
GJ3 3 


Figure 3.65 X BAR CHART vs. CORRELATION COEFFICIENT 


AND SAMPLE SIZE 






3.12.4 UNIVARIATE X CUSUM CHART : 


From the design of experiment it can be observed that following factors are significant 

• Sample Size 

• Deviation 


Table 3.32 Summary of Effects of X_Cusum Chart 


Effect 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 

54 42222 

144 

49.90556 

1 090504 

0 338806 

2 

2 

319 0722 

144 

49.90556 

6.393521 

0 002187 

3 

3 

471 6593 

144 

49 90556 

9 451037 

9 68E-06 

12 

4 

58 44722 

144 

49 90556 

1.171157 

0 326057 

13 

6 

54 1037 

144 

49 90556 

1 084122 

0 374739 

23 

6 

51 97593 

144 

49.90556 

1.041486 

0 401012 

123 

12 

46 66204 

144 

49 90556 

0 935007 

0.514044 


Following are the observations made from the results as shown in Figure 3.66 to 3 70 

• With the increase in correlation coefficient, this chart performance declines but 
after second level, performance again starts to improve. 

• As sample size n increases, performance increases. But after sample size n = 6, 
improvement in performance is slow. 

• At high deviation, signals are faster than low deviation. 
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Variable X_CUSUM 4X' Vanable X CUSUM 


Plot of Means (unweighted) 
SAM_S1ZE Mam Effect 
F(2,144)=6 39, p< 0022 



Plot of Means (unweighted) 
DEVI ATI O Mam Effect 
F(3,144) s =9 45, p< 0000 



;ure 3.67 X_CUSUM vs. DAVIATION 


Plot of Means (unweighted) 
COR_COF Mam Effect 
F(2,144)-1 09, p< 3388 



Figure 3.68 X_CUSUM vs. CORRLATION COEFFICIENT 






3.13 Design of Experiment for Five Variables with Run-Length as Performance 
Measure 

Performance of charts is also examined on five variables to see the impact of deviation 
w.r.t to number of variable. This design of experiment is also conducted for 3 levels of 
correlation coefficient, 3 levels of sample size and 4 levels for deviations. The run length 
is taken as output for each combination of these factors. 

Following are the factors and levels for this design of experiment. 


Table 3.33 Factors of DOE 


Factor 

Levels 

Correlation Coefficient 

Low, Middle and High 

Sample Size 

1,6 and 10 

Deviation 

Low, Middle, High and Combination 

of Low and High 


In case of five variables, the correlation coefficient between variables are discussed in the 
section-3.1. The levels of deviation are alike two variable. 

3.14 Analysis of DOE 

Six charts are analyzed for 5 variable case. This experiment is replicated 7 times for all 
charts. Thus total of 256 (3 * 3 * 4 * 7) problems are solved for these charts as described 
below. 

3.14.1 T 2 CONTROL CHART : 

From the design of experiment it can be observed that following factors are significant: 

• Correlation Coefficient 

• Sample Size 


• Deviation 








• Correlation Coefficient and Sample Size 

• Correlation Coefficient and Deviation 

• Sample Size and Deviation 

• Correlation Coefficient, Sample Size and Deviation 


Table 3.34 


Effect 


Summary of Effects of T 2 Control Chart 



df 

Effect 


2 


2 


MS 

Effect 


134.9184 


263.0822 


29.8489 


21.89553 


18.07892 


216 

6.775215 

4.405602 

0.001917 

216 

6.775215 

3.23171 

0.004628 

216 

6.775215 

2.668391 

0.016154 

216 

6.775215 

1.891261 

0.036709 



































Plot of Means (unweighted) 
SAM.SIZE Mam Effect 
F(2,21 6)=38 83, p< 0000 



Figure 3.70 T 2 CHART vs. SAMPLE SIZE 


Plot of Means (unweighted) 
DEVIATIO Mam Effect 
F(3,216)=57 73, p< 0000 



Figure 3.71 


T 2 CHART vs. DEVIATION 





Variable T2_CHART 


Plot of Means (unweighted) 
2-way interaction 
F(6,216)=2 67, p< 0162 



-o- SAM.SIZE 
G_1 1 

□ SAM SIZE 
G_2 2 

«■ SAM.SIZE 
G_3 3 


Figure 3.72 T 2 CHART vs. DEVIATION AND SAMPLE SIZE 


Plot of Means (unweighted) 
2-way interaction 
F(6,216)=3 23, p< 0046 



~ o — COFLCOFF 
GJ 1 

□ COFLCOFF 
G_2 2 

❖ COR^COFF 
G_3 3 


Figure 3.73 T 2 CHART vs. DEVIATION AND CORRELATION 

COEFFICIENT 





3 way interaction 
F{ 1 2,21 6)=1 89, p< 0367 



—o— COFLCOFF 
G_1 t 

□ COR_COFF 
G_2 2 

o COR_COFF 
G_3 3 


Figure 3.74 T 2 CHART vs. DEVIATION, CORRELATION 


COEFFICIENT AND SAMPLE SIZE 


3.14.2 MULTIVARIATE CUSUM CHARTS 

3.14.2.1 CUSUM #1: 

From the design of experiment it can be observed that following factors are significant: 

• Correlation Coefficient 

• Sample Size 

• Deviation 

• Correlation Coefficient and Sample Size 


Table 3.35 Summary of Effects for Cusum #1 


Effects 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 

95.91521 

216 

4.576279 

20.95921 

4.8E-09 

2 

2 


216 

4.576279 

13.04124 

| 

3 

3 

257.731 1 

216 

4.576279 

56.31892 

6.18E-27 

12 

4 

30.72733 

216 

4.576279 

6.71448 

4.12E-05 

13 

6 

5.515128 

216 

4.576279 

1.205156 

0.304698 

23 

6 


216 

4.576279 

0.407892 

0.873355 

123 

12 

2.555474 

216 

4.576279 

0.558417 

0.873637 










































GJ'I G_2 2 G_3,3 

SAM_SIZE 


Figure 3.76 CUSUM#1 CHART vs. SAMPLE SIZE 




Plot of Means (unweighted) 
DEVIATIO Main Effect 
F(3,21 6)=56 32, p< 0000 



Figure 3.77 CUSUM#1 CHART vs. DEVIATION 


Plot of Means (unweighted) 
2-way interaction 
F(4,216)=6 71, p< 0000 



-<►- COR_COFF 
G_1 1 

□ COR_COFF 
G_2 2 

COR.COFF 
G_3 3 


Figure 3.78 CUSUM #1 CHART vs. SAMPLE SIZE AND 


CORRELATION COEFFICIENT 


3.14.2.2 CUSUM #2: 


From the design of experiment it can be observed that following factors are significant: 


Correlation Coefficient 
Sample Size 





• Deviation 


• Correlation Coefficient and Sample size 

• Correlation Coefficient and Deviation 


Table 3.36 Summary of Effects for Cusum#2 


Effect 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 


216 

2.61993 

85.48335 

4.49E-28 

2 

2 

29.95492 

216 

2 61993 

11.43348 

1.91E-05 

3 

3 

156.8853 

216 

2.61993 



12 

4 

7.288027 

216 

2.61993 

2.781765 

0.027713 

13 

6 


216 

2.61993 

16.39251 

1.53E-15 

23 

6 

3.444307 

216 

2.61993 

1.314656 

0.251695 

123 

12 

4.69612 

216 

2.61993 

1.79246 

0.050775 




























Plot of Means (unweighted) 
COR.COFF Main Effect 
F(2,216)=85 48, p< 0000 



Figure 3.79 CUSUM #2 vs. CORRELATION COEFFICINT 


Plot of Means (unweighted) 
SAM.SIZE Mam Effect 
F(2,21 6)=1 1 43, p< 0000 



Figure 3.80 


CUSUM #2 vs. SAMPLE SIZE 




Variable MC2_CHAR 








Plot of Means (unweighted) 
2-way interaction 
F(6,216)=t6 39 , p< 0000 



Figure 3.83 CUSUM#2 vs. CORRELATION 
COEFFICEINT AND DEVIATION 


— o— COR_COFF 
G_1 1 

□ COR_COFF 
G_2 2 

« COR_COFF 
G_3 3 


3.14.3 Z.METHOD : 


From the design of experiment it can be observed that following factors are significant: 

• Correlation Coefficient 

• Deviation 

• Correlation Coefficient and Sample Size 

• Correlation Coefficient and Deviation 


Table 3.37 Summary of Effects for Z_Method 


Effect 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 


216 



0.00069 

2 

2 

7.045756 

216 



0.325105 

3 

3 

173.0808 

216 

6.238095 



12 

4 

35.16052 

216 

6.238095 

5.636419 

0.000247 

13 

6 

18.18348 

216 




23 

6 

1.665526 

216 

6.238095 

0.266993 

0.9518 

123 

12 

3.417421 

216 

6.238095 

0 547831 

0.881389 






























Variable Z_METHOD 


Plot of Means (unweighted) 
COR.COFF Main Effect 
F(2,216)=7 53, p< 0007 



Figure 3.84 ZJMETHOD vs. CORRELATION COEFFICIENT 


Plot of Means (unweighted) 
SAM_SIZE Mam Effect 
F(2,216)-1 13, p< 3251 



Figure 3.85 


Z_METHOD vs. SAMPLE SIZE 




Variable Z_METHOD 


Plot of Means (unweighted) 
DEVIATIO Mam Effect 
F(3,216)=27 75, p< 0000 



DEVIATIO 


Figure 3.86 


ZJMETHOD vs. DEVIATION 


Plot of Means (unweighted) 
2-way interaction 
F(4,21 6)-5 64, p< 0002 



—o— COR,COFF 
G J 1 

□ COR,COFF 
G_2 2 

o COR,COFF 
G,3*3 


SAM, SIZE 


Figure 3.87 Z_METHOD vs. SAMPLE SIZE AND 

CORRELATION COEFFICEINT 




3.14.3 SHEWHART X_BAR CHART 


From the design of experiment it can be observed that following factors are significant: 

• Correlation Coefficient 

• Deviation 

• Correlation Coefficient and Sample Size 

• Correlation Coefficient and Deviation 

Table 3.38 Summary of Effects for XJBar Chart 


Effect 

df 

Effect 

MS 

Effect 


m 

F 

p-level 

1 

2 

9.133263 

216 

6.308422 

1 .447789 

0.237362 

2 

2 

347.6745 

216 

6.308422 

55.11276 

4.58E-20 

3 

3 

517.1519 

216 

6.308422 

81.97803 

0 

12 

4 


216 

6.308422 

1.033487 

0.390828 

13 

6 

2.358357 

216 

6.308422 

0.373843 


23 

6 

14.8939 

216 

6.308422 

2.360956 

0.031398 

■a 

12 

j|j|g|||j 

216 





Plot of Means (unweighted) 
COR.COFF Main Effect 
F(2,216)s=1 45, p< 2374 



Figure 3.88 X_BAR CHART vs. CORRELATION COEFFICIENT 

















Variable X BAR 







3.14.4 UNIVARIATE X_CUSUM CHART : 


From the design of experiment it can be observed that following factors are significant: 

• Sample Size 

• Deviation 


Table 3.39 Summary of Effects for X_Cusum Chart 


Effect 

df 

Effect 

MS 

Effect 

df 

Error 

MS 

Error 

F 

p-level 

1 

2 

3.383108 

216 

5.118937 

0.6609 

0.517427 

2 

2 

220.1268 

216 

5.118937 

43.00245 

1.9E-16 

3 

3 


216 



3.7E-28 

12 

4 

2.328756 

216 

5.118937 

0.45493 

0.768737 

13 

6 

1 .083936 

216 

5.118937 

6.21175 


23 

6 

4.596652 

216 

5.118937 

0.89797 

0.497215 

123 

12 

1 .626521 

216 

5.118937 

0.317746 






















Figure 3.93 X_CUSUM CHART vs. DEVIATION 






CHAPTER 4 


RESULTS AND CONCLUSIONS 

This chapter discusses the results obtained after the complete analysis of all charts for a 
number of problems and the conclusions that can be drawn form the analysis. 

4.1 Results and Conclusions 

• Charts which is basically for detecting process mean deviation or process variance 
deviation gives signal when process variance deviation or process mean deviation is 
introduced 

• Performance of following charts increases as correlation coefficient increases 

1 T 2 Control Chart 

2 x 1 Control Chart 

3 Cusum #1 Chart 

4 Z_Method 

5 X_Bar Chart 

6 Exponentail_Chart 

• Performance of Cusum #2 Chart initially increases as correlation coefficient 
increases but aftei 0.7 p , its becomes constant 

• Peiformance of unr variate X_Cusum chart declines as correlation coefficient 
increases 

• Peiformance of following chart improves as sample size increases 

1. Cusum #1 Chart 

2. X_Bar Chart 

3 Univariate X_Cusum 


123 



® Performance of following chart improves as sample size increases but after sample 
size n = 6, its becomes constant 

1 T 2 Control Chart 

2 x 1 Control Chart 

3 Cusum#2 Chart 

• Pei formance of following chart declines as sample size increases 

1 Z_Method 

2 Exponential_Chart 

• For every chart, as deviation increases, they perform better 

• All charts perform better in positive process mean or process variance deviation than 
same value of negative process mean or process variance deviation 


4.2 Scope of the Future Work 

In the the present work an attempt is made to study the performance of different 
multivariate quality-control charts w r t correlation coefficient, sample size and process 
mean or variance deviation. 

Although these charts have triggering capabilities that process is going out-of-control but 
only T 2 chart and Z_method have diagnostic procedure for out-of-control variables 
Others do not have diagnostic procedure, which can be a big problem if numbers of 
variables are large Multivariate Cusum Charts give good results but its mam drawback is 
the absence of diagnostic procedure for out-of-control signal A study may be done for 
developing the diagnostic procedure of these charts. 
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